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I.  INTRODUCTION: 

Previous  MRI  studies  in  multiple  sclerosis  (MS)  have  retrospectively  analyzed  normal-appearing 
brain  tissue  in  locations  where  typical  MS  lesions  ultimately  appeared,  finding  pre-lesional 
changes  in  several  MRI  metrics.  However,  studies  have  not  been  entirely  consistent  and  the 
development  of  a  prototypical  MS  lesion  cannot  as  yet  be  prospectively  predicted.  The  primary 
objective  of  this  project  is  to  validate  the  “preactive”  lesion  hypothesis  in  MS  by  identifying  the 
spatiotemporal  imaging  signature  of  white  matter  (WM)  destined  to  undergo  acute,  focal 
inflammation  and  demyelination — specifically,  one  that  will  allow  reliable,  prospective  detection 
of  nascent  lesions  before  they  appear  on  conventional  (non-quantitative)  imaging.  The  specific 
aim  is  to  acquire  a  longitudinal  set  of  quantitative  MRI  metrics  in  MS  patients  and  perform  a 
multivariate  spatiotemporal  analysis  of  pre-lesional,  normal-appearing  white  matter,  seeking 
spatially  clustered  interval  changes  that  presage  the  appearance  of  a  typical  MS  plaque. 

II.  BODY: 

The  key  accomplishments  achieved  over  the  course  of  the  project  are  detailed  below,  keyed  to 
the  tasks  outlined  in  the  Statement  of  Work. 


II.A.  Subject  enrollment: 

Recruitment  for  the  study  was  challenging  as  several  disease-modifying  therapies  were  approved 
by  the  FDA  and  made  available  to  patients  during  this  phase  of  the  study;  this  made  it  difficult  to 
find  subjects  with  active  disease  (necessary  for  new  lesions  to  study)  who  were  willing  to  forego 
therapy  for  the  duration  of  the  study  (required  by  inclusion  criteria).  Seven  treatment-naive  MS 
patients  meeting  inclusion  criteria  were  ultimately  enrolled  and  completed  the  eight  monthly 
MRI  scans.  As  it  happened,  the  delay  in  enrollment  may  have  been  a  “blessing  in  disguise”  as  it 
allowed  us  to  continue  the  preparatory  development  and  optimization  of  MRI  acquisition 
methods,  leading  to  several  key  accomplishments  as  detailed  in  the  next  section. 


II.B.  MRI  acquisition: 

We  made  several  technical  advances  in  quantitative  MRI  acquisition  methods  for  the  project, 
specifically  in  quantitative  magnetization  transfer  (MT)  imaging  [1-7]. 


We  developed  a  novel  method  for  simultaneous  calibration  of  excitation  flip  angle  and  MT 
saturation  power  [1-2];  accounting  for  B1  and  flip  angle  variations  is  critical  for  longitudinal 
stability  of  MT  imaging  which,  in  turn,  is  required  for  the  tracking  of  changes  that  occur  in 
lesions  over  time.  Previous  studies  relied  on  either  B1  maps  or  flip  angle  maps,  but  not  both,  to 
correct  both  excitation  flip  and  MT  saturation  power,  which  may  result  in  significant  errors  in 
MT  parameters.  Our  new  method  enabled  the  accounting  for  both  sources  of  errors  and,  as  a 
result,  yielded  more  accurate  MT  measures  (Fig.  1). 


Figure  1.  Left:  Bound  pool  fraction  /  in  a 
homogeneous  gel  phantom  calculated  using  old  MT 
method  AFI  (FA)  and  new  method  (BS-AFI).  Right: 
Representative  profiles  through  the  images  showing  the 
old  method  overestimating  the  key  MT  measure  in  the 
areas  close  to  the  edges. 
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We  also  developed  a  new  quantitative  MT  imaging  method,  modified  cross-relaxation  imaging, 
(mCRI)  [5],  which  decreases  the  number  of  acquired  images  and  scan  time  relative  to  existing 
methods  while  yielding  more  accurate  (by  -16%)  values  of  macromolecular  proton  fraction 
( MPF ),  an  imaging  biomarker  of  myelin  content.  This  was  achieved  via  novel  modeling  of  MT 
exchange  processes  in  steady-state  sequences,  an  existing  method  that  was  transformed  into  a 
fast,  minimum-acquisition  protocol  capable  of  mapping  a  subset  of  qMTI  measures  including 
MPF  and  T1  using  only  three  MRI  scans  [6-7]  (up  to  60%  reduction  over  state-of-the-art 
techniques).  This  advance  was  achieved  by  applying  physiologically  relevant  constraints  on  the 
remaining  model  parameters  and  optimizing  acquisition  to  minimize  the  impact  of  these 
parameters  on  accuracy.  Fig.  2  shows  results  of  fast  mCRI  and  illustrates  the  limitations  of 
R1=\/T1  as  an  indirect  measure  of  myelin.  Fig.  3  favorably  compares  MPF  with  an  alternative 
MT-based  myelin  marker  based  on  inhomogeneous  MT  effect  ( ihMTR ).  The  final  7  min  version 
of  the  protocol  was  applied  in  study  subjects. 

4 
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Figure  2.  MPF  and  R1=1/T1  maps  using  our  fast 
three-measurement  mCRI  with  flip  angle  correction. 
Note  excellent  myelin  contrast  on  MPF  map. 


MPF  imaging  in  WM  is  illustrated  by  Fig.  5  below.  Note  that  two  different  lesions  appearing 
similar  on  standard  T2-FLAIR  images  demonstrate  different  levels  of  tissue  damage  on  MPF 
(green  arrow  -  lesion  with  high  residual  myelination,  red  arrow  -  lesion  with  virtually  no 
residual  myelination).  Fig.  6  illustrates  the  promise  of  MPF  to  improve  detection  of  focal  GM 
pathology  compared  to  conventional  MRI  and  state-of-the-art  double  inversion  recovery  (DIR) 


Figure  3.  Comparison  of  MPF  and  ihMTR  myelin- 
sensitive  maps  optimized  to  run  in  the  same  time  (7 
min;  ihMTR:  1.9x1. 9x3 mm,  60  slices;  MPF:  1.5  mm 
isotropic,  120  slices).  Note  much  higher 
resolution/quality  of  MPF. 


T2  FLAIR  MPF  T2  FLAIR  MPF 


Figure  5.  Two  MS  lesions  appear  similar  on 
standard  T2  FLAIR  but  different  on  MPF, 
presumably  reflecting  different  levels  of 
demyelination. 


MPF  Double  IR  IR-SPGR  T2  FLAIR 


Figure  6.  GM  lesions  on  MPF,  DIR  (state-of-the-art  GM  imaging 
method  in  MS),  and  standard  IR-SPGR  and  3D  T2  FLAIR.  Note 
improved  GM  lesion  depiction  on  MPF. 
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II.C.  Image  Processing: 

We  developed  a  semi-automated  image  processing  pipeline  for  multiparametric  images,  with  key 
stages  including  brain  extraction,  co-registration  of  images  from  different  modalities  and  time- 
points,  and  segmentation  of  normal-appearing  white  matter  (Figs.  7-8).  All  imaging  data 
acquired  in  different  orientations  (axial  and  sagittal)  were  first  converted  to  the  NIFTI  image 
format  and  their  orientations  changed  to  axial  (orientation  of  the  reference  data).  Skull  stripping 
(brain  extraction)  was  then  performed  on  high-resolution  3D  BRAVO  T1 -weighted  images  using 
the  Brain  Extraction  Tool  in  the  FSL  software  library  (FMRIB,  Oxford,  UK).  To  correct 
intensity  variations  due  to  magnetic  field  inhomogeneity,  non-uniformity  correction  was  applied 
on  the  extracted  brain  BRAVO  image  using  JIM  analysis  software  (Xinapse  Systems,  UK).  This 
reference  BRAVO  image  was  then  co-registered  to  the  MNI152  template  (Montreal 
Neurological  Institute)  using  the  FLIRT  (FMRIB  Linear  Image  Registration  Tool)  brain  image 
registration  toolbox  in  FSL.  Eddy  current  correction  was  applied  to  DTI  data  also  using  FSL 
software.  All  images  for  a  given  subject  at  all  scan  dates  were  co-registered  to  that  subject’s 
reference  BRAVO  series  and  the  BRAVO  brain  mask  generated  in  the  first  step  was  applied  to 
all  co-registered  images.  In  addition,  non-uniformity  correction  was  applied  to  high-resolution 
3D  Cube  T2-FLAIR  images. 
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Figure  7:  Overview  of  post-processing  pipeline 
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Figure  8:  Block  diagram  of  image  co-registration  steps  to  align  quantitative  (DTI,  MT)  with  anatomical  (Tl,  T2)  images. 


II.D.  Data  Analysis: 

The  image  processing  detailed  above  yielded  a  set  of  spatially  and  temporally  co-registered, 
lesion-segmented,  quantitative  MRI  maps  that  could  be  analyzed  to  identify  any  pre-lesional 
changes  in  MR  parameters.  An  example  is  shown  in  Fig.  9.  Lesions  that  appeared  for  the  first 
time  in  any  of  the  second-half  MRI  scans  (i.e.  with  >4  pre-lesional  scans  available)  were  selected 
for  analysis,  seeking  statistically  significant  differences  in  MRI  parameters  between  the  NAWM 
that  evolved  into  a  lesion  versus  the  NAWM  that  did  not. 


Gd+TI 


Figure  9.  Serial  monthly  magnetization  transfer  ratio  (MTR)  maps  spatiotemporally  co-registered  to  post-contrast  Tl- 
weighted  images  (Gd+TI)  for  a  subject  who  developed  a  new  MS  lesion  in  month  #8  (arrows).  Two  regions  of  interest 
(ROIs),  one  encompassing  the  lesion  and  another  normal-appearing  white  matter  (NAWM),  could  then  be  tracked  over  the 
7  monthly  scans  before  the  lesion  appeared.  Note  that  MTR  in  the  lesion  (blue)  decreased  precipitously  relative  to  NAWM 
(red)  when  the  lesion  appeared  but  not  before  then  (ratio  of  lesion  to  NAWM  in  green). 


Analyzed  parameters  included  four  HYDI  metrics:  fractional  anisotropy  (FA),  axial  (Da)  and 
radial  (Dr)  diffusivity  and  zero-displacement  probability  (Po),  and  four  MT  metrics: 
magnetization  transfer  ratio  (MTR),  MT  saturation  (MTSat),  bound  pool  fraction  (f)  and 
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longitudinal  relaxation  rate  (Rl).  There  were  no  statistically  significant  differences  between  pre- 
lesional  WM  and  other  NAWM  on  any  of  the  pre-lesional  scans,  for  any  of  the  parameters 
analyzed,  that  would  have  allowed  the  prediction  of  lesion  appearance  beforehand. 


Fractional  Anisotropy  (FA)  Axial  Diffusivity  (Da)  Bound  Pool  Fraction  (f)  Longitudinal  Relax.  Rate  (Rl) 


1.80  E-03 

1.40E-01 

1.20E01 

1.40  E*- 00 

t  1.20  E+ 00 

-  1.00  E+ 00 

hm  j  T - ^ - 1 

1  1  \  1.20  E-03 

\  1.00  E-03 

V  8.00  E  04 

\  6.00E-04 

1  i  l  ^ 

I  8.00  E-02 

6.00E-02 

8.00  E-01 

- -  6.00  E-01 

—  A 

4.00E-04 

4.00E-02 

2.00  E-02 

- i —  4.00  E-01 

2.00  E-01 

„  0.00  E+ 00 

4-3-2  1  0 

SCAN 

0.00E+00 

4  3  2  -1  0 

SCAN 

0.00  E+  00 

4  -3  -2  -1  0 

SCAN 

4  -3  -2  -1  0 

SCAN 

Radial  Diffusivity  (Dr) 


Zero  Displacement  (PO) 


Mag.  Transfer  Ratio  (MTR)  Mag.Transfer  Saturation  (MTSat) 


Figure  10.  Summary  of  prelesional  diffusion  (HYDI)  and  magnetization  transfer  (MT)  parameter  tracking.  Lesion  appearance  is  on  “Scan  0.” 


Parameters  at  the  lesion  location  (blue)  are  tracked  over  the  4  previous  scans  along  with  regional  NAWM  that  does  not  evolve  into  a  lesion  (red). 
Even  on  the  scan  immediately  preceding  lesion  appearance  (“Scan  -1”)  there  were  no  significant  changes  in  any  parameter  that  would  allow 
lesion  development  to  be  anticipated  beforehand. 


III.  KEY  RESEARCH  ACCOMPLISHMENTS: 

•  Developed  novel  method  for  simultaneous  calibration  of  excitation  flip  angle  and  MT 
saturation  power,  thus  improving  the  accuracy  and  longitudinal  stability  of  MT  imaging 
[1-2]. 

•  Developed  novel  quantitative  MT  imaging  method,  “modified  cross-relaxation  imaging” 
(mCRI)  [5],  which  decreases  the  number  of  acquired  images  and  thus  scan  time  relative 
to  existing  methods  while  yielding  more  accurate  (by  ~16%)  values  of  macromolecular 
proton  fraction,  an  imaging  biomarker  for  myelin. 

•  Developed  semi-automated  image  processing  pipeline  for  longitudinal  studies  of 
multiparametric  images,  with  key  stages  including  co-registration  across  different 
modalities  and  time-points  as  well  as  tissue  segmentation.  This  pipeline  will  be  highly 
useful  for  any  future  studies  involving  longitudinal  tracking  of  multiparametric  imaging 
data. 

•  Achieved  monthly  tracking  of  multiple  MRI  metrics  in  treatment-naive  MS  patients  with 
active  disease,  allowing  quantitative  analysis  of  the  natural  history  of  lesion  development. 

IV.  REPORTABLE  OUTCOMES: 

IV. A.  Manuscripts: 

•  Hurley  SA,  Yamykh  VL,  Johnson  KM,  Field  AS,  Alexander  AL,  Samsonov  AA. 
Simultaneous  variable  flip  angle-actual  flip  angle  imaging  method  for  improved  accuracy 
and  precision  of  three-dimensional  T1  and  B1  measurements.  Magn  Reson  Med 
2012;68(l):54-64. 
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•  Mossahebi  P,  Yamykh  VL,  Samsonov  A.  Analysis  and  correction  of  biases  in  cross¬ 
relaxation  MRI  due  to  biexponential  longitudinal  relaxation.  Magn  Reson  Med 
2014;71(2):830-838. 

IV.B.  Abstracts: 

•  Hurley  SA,  Mossahebi  P,  Johnson  KM,  Samsonov  AA.  Simultaneous  Mapping  of  B1 
and  Flip  Angle  by  Combined  Bloch-Siegert,  Actual  Flip-Angle  Imaging  (BS-AFI).  In: 
Proceedings  of  the  International  Society  for  Magnetic  Resonance  in  Medicine  (ISMRM) 
20th  Scientific  Meeting,  Melbourne,  Australia,  May  2012. 

•  Mossahebi  P,  Yarknykh  VL,  Samsonov  AA.  Improved  Accuracy  of  Cross-Relaxation 
Imaging  Using  On-Resonance  MT  Effect  Correction.  In:  Proceedings  of  the  International 
Society  for  Magnetic  Resonance  in  Medicine  (ISMRM)  20th  Scientific  Meeting, 
Melbourne,  Australia,  May  2012. 

•  Mossahebi  P,  Samsonov  AA.  Optimization  Strategies  for  Accurate  Quantitative  MT 
Imaging.  In:  Proceedings  of  the  International  Society  for  Magnetic  Resonance  in 
Medicine  (ISMRM)  20th  Scientific  Meeting,  Melbourne,  Australia,  May  2012. 

•  Samsonov  AA,  Mossahebi  P,  Anderson  A,  Velikina  JV,  Johnson  KM,  Johnson  SC, 
Fleming  JO,  Field  AS.  High  Resolution,  Motion  Corrected  Mapping  of  Macromolecular 
Proton  Fraction  (MPF)  In  Clinically  Acceptable  Time  Using  3D  Undersampled  Radials. 
In:  Proc  of  ISMRM  2014;  Milan,  Italy,  p  3337. 

•  Mossahebi  P,  Samsonov  A.  Rapid  and  Accurate  Variable  Flip  Angle  T1  Mapping  with 
Correction  of  MT  Effects.  In:  Proceedings  of  the  International  Society  for  Magnetic 
Resonance  in  Medicine  (ISMRM)  20th  Scientific  Meeting,  Melbourne,  Australia,  May 
2012. 

IV.C.  Patents: 

•  “Efficient  and  Deterministic  Approach  to  Generating  Optimal  Ordering  of  Diffusion  MRI 
Measurements.”  Patent  applied  for  and  awarded  to  Samuel  A.  Hurley,  the  graduate 
student  who  worked  on  the  project  and  was  supported  by  the  award  during  his 
development  of  the  invention.  This  is  a  novel  ordering  scheme  for  diffusion  encoding 
gradients  used  to  acquire  diffusion-weighted  MR  images.  The  invention  achieves 
increased  efficiency  of  image  acquisition  several  orders  of  magnitude  greater  than 
previously  existing  schemes. 

IV. D.  Degrees  and  Employment  Opportunities 

The  graduate  student  who  worked  on  the  project,  Samuel  A.  Hurley,  received  his  Ph.D.  in 
Medical  Physics  while  supported  by  the  award.  He  also  successfully  applied  for  a  post-doctoral 
research  scientist  position  at  the  University  of  Oxford  Centre  for  Functional  Magnetic  Resonance 
Imaging  of  the  Brain  (FMRIB). 

V.  CONCLUSION: 

The  specific  aim  of  the  project  was  to  acquire  a  longitudinal  set  of  quantitative  MRI  metrics  in 
MS  patients  and  perform  a  multivariate  spatiotemporal  analysis  of  pre-lesional,  normal¬ 
appearing  white  matter,  seeking  spatially  clustered  interval  changes  that  presage  the  appearance 
of  a  typical  MS  plaque.  Some  previous  studies  have  found  such  changes  while  others  have  not, 
and  to  date  there  is  no  known  method  for  reliably  predicting  the  location  of  a  plaque  before  it  is 
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visibly  manifest.  We  did  not  find  any  pre-lesional  changes  that  would  make  such  prediction 
possible.  Assuming  such  changes  are  indeed  occurring,  there  are  several  possible  reasons  for 
failing  to  detect  them  which  might  be  addressed  in  future  work:  (1)  changes  occurring  on  a 
spatial  scale  below  current  limits  of  MRI  resolution;  (2)  changes  occurring  on  a  time  scale 
shorter  than  the  monthly  intervals  used  for  this  study;  or  (3)  changes  occurring  of  a  magnitude 
not  exceeding  technical  inter-scan  variance  in  MRI  measurements,  all  of  which  are  inherently 
“noisy”  to  some  degree  even  with  state-of-the-art  acquisition  and  processing  methods. 

Although  we  were  unable  to  achieve  a  prediction  of  MS  lesions  before  they  are  visibly  manifest 
(which  would  have  been  a  major,  paradigm-shifting  advance  in  the  management  of  MS  patients), 
we  did  achieve  several  significant  improvements  in  quantitative  MRI  methods  important  to  the 
study  of  MS  that  will  benefit  future  work  in  the  field. 
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Simultaneous  Variable  Flip  Angle-Actual  Flip  Angle 
Imaging  Method  for  Improved  Accuracy  and  Precision 
of  Three-Dimensional  T1  and  B1  Measurements 

Samuel  A.  Hurley,1  Vasily  L.  Yarnykh,2  Kevin  M.  Johnson,1  Aaron  S.  Field,3 
Andrew  L.  Alexander,1,4  and  Alexey  A.  Samsonov3* * 


A  new  time-efficient  and  accurate  technique  for  simultaneous 
mapping  of  T-,  and  B-,  is  proposed  based  on  a  combination  of 
the  actual  flip  angle  (FA)  imaging  and  variable  FA  methods.  Vari¬ 
able  FA-actual  FA  imaging  utilizes  a  single  actual  FA  imaging 
and  one  or  more  spoiled  gradient-echo  acquisitions  with  a  si¬ 
multaneous  nonlinear  fitting  procedure  to  yield  accurate  T1/B1 
maps.  The  advantage  of  variable  FA-actual  FA  imaging  is  high 
accuracy  at  either  short  T ^  times  or  long  repetition  times  in  the 
actual  FA  imaging  sequence.  Simulations  show  this  method  is 
accurate  to  0.03%  in  FA  and  0.07%  in  T ^  for  ratios  of  repetition 
time  to  T1  time  over  the  range  of  0.01-0.45.  We  show  for  the 
case  of  brain  imaging  that  it  is  sufficient  to  use  only  one  small 
FA  spoiled  gradient-echo  acquisition,  which  results  in  reduced 
spoiling  requirements  and  a  significant  scan  time  reduction 
compared  to  the  original  variable  FA  method.  In  vivo  validation 
yielded  high-quality  3D  7^  maps  and  T ^  measurements  within 
10%  of  previously  published  values  and  within  a  clinically  ac¬ 
ceptable  scan  time.  The  variable  FA-actual  FA  imaging  method 
will  increase  the  accuracy  and  clinical  feasibility  of  many  quan¬ 
titative  MRI  methods  requiring  T^/B^  mapping  such  as  dynamic 
contrast  enhanced  perfusion  and  quantitative  magnetization 
transfer  imaging.  Magn  Reson  Med  68:54-64,  2012.  ©  2011 
Wiley  Periodicals,  Inc. 

Key  words:  T,  relaxometry;  B A  mapping;  actual  flip  angle 
imaging  (AFI);  spoiled  gradient  echo 


Accurate  measurements  of  longitudinal  relaxation  time 
(7\)  are  essential  to  many  quantitative  MRI  techniques 
and  clinical  applications.  Traditional  inversion  recovery 
T1  mapping  methods  often  result  in  prohibitively  long 
acquisitions  due  to  a  long  repetition  time  (TR)  (1).  The 
Look-Locker  method  (2)  may  be  used  to  accelerate  the 
inversion  recovery  experiment  by  frequent  sampling  of 
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the  entire  recovery  curve.  This  two-dimensional  (2D) 
approach  can  be  modified  to  allow  multislice  volumetric 
coverage  (3,4).  A  three-dimensional  (3D)  Look-Locker 
method  (5)  has  also  been  proposed  to  improve  the  sig- 
nal-to-noise  ratio  (SNR)  and  slice  direction  resolution  of 
2D  techniques  at  the  cost  of  blurring  in  phase  encode 
directions  due  to  segmented  k-space  acquisition. 

The  variable  flip  angle  (VFA)  method  (6,7),  also  known 
as  driven-equilibrium  single-pulse  observation  of  T1  (DES- 
POTl)  (7),  uses  several  short  TR  spoiled  gradient-echo 
(SPGR)  acquisitions  with  varying  flip  angle  (FA)  to  mea¬ 
sure  7V  It  has  gained  significant  popularity  over  the  past 
decades  due  to  its  superior  time  efficiency,  allowing  rapid 
and  accurate  3D  T1  mapping  with  high  resolution  and 
large  spatial  coverage  (8-10).  Although  the  acceleration  of 
any  T1  mapping  method  is  possible  by  the  use  of  fast 
readout  sequences  (e.g.,  echo-planar  or  spiral)  (11,12), 
applications  of  such  techniques  are  limited  to  a  large 
extent  due  to  image  artifacts  inherent  to  long  gradient- 
echo  signal  readouts.  In  contrast,  VFA  is  the  only  fast  T1 
mapping  method  to  date  that  affords  clinically  reasonable 
scan  times  and  3D  spatial  coverage  without  a  multiecho 
readout  and,  therefore,  provides  an  optimal  solution  for 
high-resolution  anatomic  applications.  In  spite  of  these 
benefits,  the  VFA  method  suffers  from  a  strong  depend¬ 
ence  on  an  accurate  knowledge  of  excitation  FA,  which 
leads  to  significant  errors  in  the  presence  of  a  nonuniform 
radiofrequency  excitation  field  (Ra),  slab  excitation  pulse 
profiles,  or  FA  miscalibration  (9,10,13).  This  is  particu¬ 
larly  problematic  at  high  field  strengths,  where  additional 
wave  effects  in  tissue  may  cause  B1  variations  up  to  30% 
at  3  T  and  in  situ  FA  calibration  and  correction  is  neces¬ 
sary  on  a  per-subject  basis  (14). 

It  is  therefore  commonly  accepted  that  VFA  T1  map¬ 
ping  at  3  T  or  higher  field  strengths  needs  to  be  com¬ 
bined  with  an  appropriate  B1  correction  method  (9,10), 
and  many  methods  to  measure  FA  have  been  developed 
and  applied  to  the  correction  of  VFA  T1  measurements. 
Similar  to  inversion  recovery-based  T1  mapping  techni¬ 
ques,  the  double  angle  method  (15)  uses  a  2D  multislice 
acquisition  with  a  long  TB  to  avoid  T1  weighting  of  the 
acquired  maps.  Scan  time  can  be  reduced  by  minimizing 
T1  dependence  with  specialized  radiofrequency  presatu¬ 
ration  pulses  (16,17)  to  achieve  a  uniform  saturation  of 
magnetization  over  a  large  area  prior  to  readout.  How¬ 
ever,  the  2D  implementation  imposes  restrictions  on  the 
use  of  double  angle  method  as  a  correction  technique  for 
3D  VFA  T1  measurements  due  to  slice  profile  effects. 
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It  is  important  to  note  that  a  3D  B1  mapping  technique 
with  matched  geometry  is  preferable  in  the  context  of 
VFA  T1  mapping  to  simplify  registration  between  data  vol¬ 
umes.  One  3D  approach  encodes  the  FA  in  the  phase  of 
the  MR  signal  (18).  This  reduces  the  7\  dependence  due 
to  previous  magnetization  history  and  can  image  over  a 
large  dynamic  range  of  FA,  but  is  sensitive  to  main  field 
inhomogeneity  and  requires  a  long  TR  relative  to  T1  to 
avoid  phase  oscillations.  Other  fast  techniques  have 
recently  been  proposed  based  on  a  stimulated-echo  echo- 
planar  acquisition  (19)  and  a  modification  of  double  angle 
method  using  a  catalyzed  preparative  sequence  (20).  Meth¬ 
ods  that  use  a  steady-state  gradient  echo  readout  are  par¬ 
ticularly  well  suited  as  these  pulse  sequences  are  closely 
related  to  the  SPGR  sequence  used  in  VFA  and  can,  there¬ 
fore,  be  acquired  with  the  same  matrix  size  and  readout 
scheme,  greatly  increasing  the  consistency  between  the  FA 
and  T1  mapping  steps.  Such  methods  include  a  multipoint 
search  of  the  signal  null  corresponding  to  the  180°  FA 
(21,22),  a  Ra-dependant  shift  in  the  spin  resonance  fre¬ 
quency  (23),  an  inversion-recovery  SPGR  termed  driven 
equilibrium  single-pulse  observation  of  T1  with  high-speed 
incorporation  of  radiofrequency  (RF)  field  inhomogeneities 
(DESPOTl-HIFI)  (24),  and  a  dual  TR  SPGR  acquisition 
termed  actual  flip  angle  imaging  (AFI)  (25). 

The  accurate  combination  of  a  R1  mapping  technique 
with  3D  VFA  T1  mapping  requires  consideration  of  addi¬ 
tional  factors  affecting  FA  values,  including  excitation 
slab/slice  profile  and  RF  pulse  properties  (26).  Two 
recently  proposed  methods,  DESPOTl-HIFI  (24)  and  AFI 
(25),  measure  FA  using  the  exact  same  excitation  pulse  en¬ 
velope  and  slab  geometry  as  the  VFA  experiment,  allowing 
measurements  with  a  very  similar  FA  distribution  present 
in  both  sequences.  DESPOTl-HIFI  requires  an  additional 
Fh-insensitive  adiabatic  inversion  preparation,  which,  with 
proper  tuning  of  the  inversion  time  to  the  expected  T1 
range,  enables  fast  and  accurate  FA  mapping.  To  avoid  Ta 
weighting,  the  AFI  method  estimates  FA  maps  using  the 
ratio  of  two  signals  from  two  different  TRs  (TRAFIlj2)  using 
a  linearized  form  of  the  signal  equation  (see  next  section), 
which  holds  under  a  main  assumption  that 

TRafii,2  <  Ti  [1] 

The  calculated  FA  map  can  then  be  applied  as  a  preca¬ 
libration  step  to  a  traditional  VFA  T1  measurement.  AFI 
is  well  suited  for  the  correction  of  VFA  measurements 
because  of  the  identical  readout  scheme  to  the  SPGR 
pulse  sequence,  the  use  of  an  identical  RF  excitation 
pulse,  and  lack  of  a  need  for  additional  magnetization 
preparation  with  specialized  RF  pulses  (20,23,24). 

In  spite  of  the  many  benefits  of  AFI,  recent  studies 
have  shown  that  large  spoiler  gradients  are  required  for 
effective  suppression  (spoiling)  of  transverse  magnetiza¬ 
tion  and  thus  accurate  FA  measurements  (27).  Such  large 
gradients  require  a  considerable  increase  in  the  TR  com¬ 
pared  to  a  standard  short  TR  SPGR  sequence  and,  there¬ 
fore,  an  overall  increase  in  scan  time.  Although  proper 
spoiling  is  essential  for  accurate  FA  quantification, 
increasing  the  TR  of  the  AFI  sequence  may  lead  to  a  vio¬ 
lation  of  the  main  assumption  (Eq.  1)  (25),  posing 
another  significant  problem  with  the  accuracy  and  over¬ 


all  efficiency  of  the  technique.  Similarly,  strong  spoiler 
gradients  are  also  necessary  in  VFA  SPGR  measurements 
for  accurate  7\  mapping,  especially  at  larger  FAs  (27,28). 
This  further  increases  the  overall  time  of  the  technique. 

An  independent  mechanism  for  violation  of  the  AFI 
assumption  (Eq.  1)  is  T1  shortening  due  to  contrast  agents, 
which  may  render  postcontrast  FA  measurements  highly 
inaccurate.  Postcontrast  T1  measurements  should  there¬ 
fore  be  corrected  with  a  precontrast  R1  map.  However, 
this  is  not  feasible  in  some  applications,  such  as  manga¬ 
nese-enhanced  MRI  of  axon  pathways,  which  requires 
injection  of  a  contrast  agent  ~  24  h  before  the  scan  (29). 
A  wide  range  of  applications  would  greatly  benefit  from  a 
rapid  and  accurate  method  to  measure  T1  and  FA  in  the 
presence  of  TVshortening  contrast  agents  (30,31). 

AFI  and  VFA  are  highly  complementary  methods;  AFI 
suffers  inaccuracies  due  to  an  assumption  about  Tlf 
whereas  VFA  suffers  inaccuracies  due  to  an  assumption 
about  FA.  In  this  article,  we  propose  an  efficient  method 
to  tackle  the  quantification  of  T1  and  R1  in  the  problem¬ 
atic  AFI  regime  (TRAFI1)2  ~  Ta)  arising  from  a  long  AFI 
TR  due  to  spoiling  requirements  and/or  short  T1  times.  In 
our  method,  variable  flip  angle-actual  flip  angle  imaging 
(VAFI),  T1  and  R1  are  mapped  simultaneously  using  an 
optimized  combination  of  AFI  and  VFA  SPGR  acquisi¬ 
tions.  By  exploiting  the  close  relationship  of  the  SPGR 
and  AFI  signal  equations,  we  demonstrate  that  the 
7^AFI1,2  ^  T\  assumption  is  not  necessary  when  T1  and 
R1  are  mapped  simultaneously,  resulting  in  a  high  FA 
mapping  accuracy  over  a  broad  range  of  T1  (32),  an 
improvement  in  the  precision  of  T1  and  R1  measurements, 
and  increased  time  efficiency  of  such  measurements.  We 
also  show  that  VAFI  reduces  the  spoiling  requirements  of 
the  SPGR  portion  of  the  method,  allowing  a  shorter  TR 
and  overall  decrease  in  acquisition  time. 

THEORY 

Variable  Flip  Angle  T -i  Mapping 

T1  measurements  based  on  VFA  SPGR  acquisitions  have 
been  previously  described  in  (6-8,33).  Given  ideal  spoil¬ 
ing  of  transverse  magnetization,  the  SPGR  signal  can  be 
described  as  follows: 

y  _ 

Sspgr  =M0  sin(a)- - — - 

1  —  cos(a) 

p  f-TR 

£i  =  exp  (it 

where  TR  and  excitation  FA  a  are  control  parameters 
prescribed  from  the  operator  console,  and  the  equilib¬ 
rium  magnetization  M0  and  T1  are  free  parameters  to  be 
determined.  Here,  we  assume  that  M0  absorbs  additional 
factors  due  to  T2*  relaxation,  receiver  gains,  and  coil 
sensitivity  profile.  In  VFA  T1  mapping,  two  or  more 
acquisitions  with  varying  a  and  fixed  TR  are  acquired 
and  fitted  to  yield  M0  and  T1.  The  problem  of  data  fitting 
can  be  conveniently  casted  into  a  linear  form  and  fitted 
with  least  squares  regression  (8)  or  the  curve  can  be  fit¬ 
ted  directly  using  nonlinear  methods,  with  the  latter 
approach  in  general  providing  better  noise  performance 
(34).  In  practice,  two  “ideal”  FAs  are  enough  to  provide 
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optimized  VFA  T1  mapping  for  a  given  T1  value  (33),  the 
lower  FA  yielding  an  image  with  proton-density  like 
contrast  and  higher  FA  with  TVlike  contrast. 

Actual  Flip  Angle  Imaging 

The  AFI  sequence  is  a  spoiled  gradient-recalled  echo  ac¬ 
quisition  with  two  identical  excitation  RF  pulses  with 
FA  a  separated  by  two  interleaved  repetitions  times 
TRafii,2  (25).  Given  complete  spoiling  of  transverse  mag¬ 
netization  at  the  end  of  each  TR  (25),  the  steady-state  sig¬ 
nal  from  the  AFI  sequence  is  described  by  the  following 
equations: 

Safii,2  =  M0sin(aAF/)x 

1  -  e[AFI2’1)  +  (1  -  ■E^n’2))£'^m’2)cos(ay4FJ) 
i  —  e[  '  ^  e[  ^cos2(aAFI)  • 
e[AFI1’  2)  =  exp(-T^m’2) 

It  was  shown  in  (25)  that  under  the  condition  of  Eq.  1 
there  exists  an  efficient  first-order  exponential  approxi¬ 
mation  without  an  explicit  dependence  on  M0  or  TV  The 
signal  then  can  be  solved  for  the  FA  as  follows: 


n  =  TRApi2  /  TRAfii  ,  r  =  SApn  /  Safii 


where  argmin(/(x))  is  the  set  of  values  x  for  which  f[x) 
attains  its  smallest  value. 

This  procedure  simultaneously  yields  a  solution  for 
three  unknown  parameters:  M0,  TV  and  k,  where  k  is  a 
scaling  factor  between  actual  and  prescribed  FA.  Theo¬ 
retically,  this  sets  the  lower  limit  on  the  total  number  of 
required  measurements  to  three.  As  a  single  AFI 
sequence  yields  two  independent  measurements,  we 
anticipate  the  minimum  number  of  SPGR  measurements 
to  be  one.  As  AFI  requires  a  large  FA  (50°-60°)  to  opti¬ 
mize  the  estimation  of  FA  maps  and  yields  images  with 
TVlike  contrast,  we  anticipate  that  the  SPGR  measure¬ 
ment  with  the  larger  FA  may  be  removed  compared  to 
the  standard  VFA  experiment  design  (8,33). 

As  the  transmit  field  varies  slowly  with  spatial  posi¬ 
tion,  an  additional  term  may  be  introduced  to  Eq.  5  to 
penalize  the  roughness  of  K(r)  and  improve  the  noise 
properties  of  reconstructed  parametric  maps  in  a  way 
similar  to  regularized  fat/water  imaging  (36),  which 
may  be  computationally  challenging  for  3D  reconstruc¬ 
tions.  To  reduce  computational  time,  we  developed  a 
simplified  three-step  version  of  the  regularization 
procedure: 

Step  1:  Solve  Eq.  5  for  spatial  maps  [M0  T\  k] 

Step  2:  Apply  spatial  smoothing  to  K(r)  to  yield  Ks(r). 
Step  3:  Solve  [M0  Ta]  =  arg  min(^Li(s^D_ 

M0  Tj 

SaFI^AFI  •  Ks))2  +  Yaa*  (sSPGR,i  ~  SsPCS{(oii.Ks))2) . 


Variable  Flip  Angle-Actual  Flip  Angle  Imaging 

To  avoid  the  main  AFI  assumption  (Eq.  1),  we  propose  a 
novel  approach,  VAFI,  which  takes  as  an  input  one  AFI 
dataset  and  one  or  more  VFA  SPGR  datasets.  In  VAFI, 
the  dependence  of  the  AFI  FA  on  T1  is  not  ignored,  but 
instead  coupled  in  a  single  nonlinear  optimization  pro¬ 
cedure  for  an  unbiased  simultaneous  estimation  of  both 
parameters.  We  assume  that  M0  is  consistent  between 
SPGR  and  AFI  acquisitions  (Eqs.  2  and  3,  respectively), 
which  requires  that  the  RF  excitation  pulse,  slab  selec¬ 
tion  gradient,  and  echo  times  be  identical  for  both  the 
SPGR  and  AFI  pulse  sequences,  and  that  the  signal  from 
both  pulse  sequences  is  perfectly  spoiled  (25,27).  Addi¬ 
tionally,  we  assume  that  error  in  a  scales  linearly  with 
transmit  field  across  the  range  of  FAs  used  in  the  VFA 
and  AFI  measurements.  This  assumption  has  been  made 
in  many  previous  FA  mapping  techniques  (15-17),  and 
recent  Bloch  equation  simulation  results  have  shown 
less  than  4%  deviation  from  linearity  over  a  range  of  2°— 
70°  (35).  For  a  set  of  SPGR  FAs  and  an  AFI  FA  aAFI 
with  corresponding  observed  signals  Sspgiu  and  SAFI1)2, 
the  unknown  parameters  can  be  determined  on  a  per- 
pixel  basis  by  simultaneously  fitting  both  the  SPGR  and 
AFI  signal  models  using  nonlinear  least  squares: 

[M0  T\  k]>  arg  min  V(sAF/,;  -  SAfiMafi  •  k)) 

M0  7\  k  \i=i 

Nspgr  \ 

+  ^2  ( SSPGR,i  -  SsPGR^i  *  k))2  J 


MATERIALS  AND  METHODS 

Implementation 

The  proposed  VAFI  fitting  method  was  implemented  in 
MATLAB  (The  Mathworks;  Natick,  MA)  using  a  built- 
in  optimization  algorithm.  As  an  initial  guess,  we  used 
the  maximum  signal  intensity  for  M0,  T1  =  Is,  and  k 
=  1.  VAFI  converged  properly  for  all  simulations  and 
experiments  without  a  need  to  modify  this  initial 
guess.  Typical  processing  time  was  ~  42  ms  per  voxel 
(64-bit  Intel  Linux  workstation,  4  cores,  32  GB  RAM). 
In  step  2  of  the  regularized  VAFI  and  standard  AFI 
approaches,  a  local  polynomial  fit  was  applied  to  gen¬ 
erate  smoothed  and  extrapolated  FA  maps  (37)  with  a 
smoothing  effect  equivalent  to  a  gaussian  kernel  with  a 
standard  deviation  of  four  voxels.  For  comparison,  T\ 
and  FA  maps  were  generated  using  standard  AFI-cor- 
rected  VFA  with  linearized  T1  estimation  and  iterative 
reweighting  (34). 

An  AFI  pulse  sequence  was  implemented  on  a  clinical 
GE  3.0T  Discovery  MR750  (GE  Healthcare;  Waukesha, 
WI)  based  on  a  product  3D  SPGR  sequence  with  the 
addition  of  a  second  TR  as  described  in  (25).  Experi¬ 
ments  were  performed  using  a  32-element  phased-array 
receive  coil  [in  vivo  experiments)  and  an  eight-element 
torso  coil  (phantom  experiments).  Spoiler  gradients  were 
designed  such  that  their  areas  AG12  were  related  as  AG2 
=  nAG1  (Eq.  4).  RF  spoiling  was  used  with  a  quadratic 
^  J  phase  increment  of  34°  (27).  For  all  experiments,  n  = 
TRafi2^TRAFi1  =  5,  as  was  suggested  in  previous  work 
(25).  The  specific  strengths  of  the  spoiler  gradients  were 
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chosen  based  on  the  completeness  of  spoiling  desired, 
herein  referred  to  as  “spoiling  regime.” 

Simulations 

Monte  Carlo  noise  simulations  were  implemented  to 
compare  the  accuracy  (bias  in  the  estimated  values)  and 
precision  (noise  in  the  estimated  values)  of  standard  two 
FA  VFA  with  AFI  correction  (henceforth  referred  to  as 
VFA)  and  VAFI.  AFI  and  SPGR  signals  were  generated 
using  Eqs.  2  and  3,  with  simulation  parameters  chosen 
to  match  the  in  vivo  experiments  described  in  the  next 
section  (SPGR:  TR  =  10  ms,  a  =  [3  18]°;  AFI:  TRAFn  = 
30  ms,  aAFI  =  55°).  Signals  were  generated  over  a  range 
of  TR/T1  =  0.01-0.45,  and  for  each  T1  time  250,000  noise 
realizations  were  created  by  adding  zero-mean  gaussian 
noise  (standard  deviation  of  le-4)  to  the  signals.  VAFI 
was  implemented  using  the  smallest  SPGR  FA  signal 
and,  to  allow  an  equivalent  comparison  with  VFA,  using 
both  SPGR  FAs.  T1  and  FA  were  then  estimated  to  assess 
the  accuracy  and  precision  of  all  three  approaches,  the 
latter  being  assessed  in  terms  of  the  T^-to-noise  and  a-to- 
noise  ratios  (7\NR  and  aNR). 

A  second  simulation  was  performed  to  investigate  the 
effect  of  FA  selection  on  the  VAFI  experiment.  For  two- 
angle  VFA,  a  pair  of  “ideal  angles"  has  been  previously 
identified  as  a  way  to  optimize  experiment  design  for  a 
specific  7\;  however,  these  angles  were  derived  based  on 
the  linearized  version  of  the  SPGR  equation  (33),  which 
is  no  longer  valid  in  the  VAFI  technique.  To  assess  the 
noise  performance  of  our  new  method,  a  set  of  45,000 
noise  realizations  (standard  deviation  of  le-4)  were  gen¬ 
erated  over  a  range  of  TR/T1  =  0.02-0.15.  In  the  first  set 
of  simulations,  aAFI  was  held  constant  at  55°  while  the 
SPGR  FA  was  varied  over  a  range  of  l°-20°.  In  the  sec¬ 
ond  set,  SPGR  FA  was  held  constant  at  5°  while  aAFI 
was  varied  over  a  range  of  40°-70°.  All  other  parameters 
match  the  previous  simulation. 

A  third  simulation  was  performed  to  evaluate  the 
effect  of  gradient  and  RF  spoiling  on  the  accuracy  of 
VAFI  quantification.  AFI  and  SPGR  signals  were  gener¬ 
ated  using  a  combined  isochromat  summation  and  diffu¬ 
sion  propagator  model  (27),  which  takes  into  account 
transverse  signal  dephasing  due  to  both  gradient  and  dif¬ 
fusion  effects.  Signals  were  generated  for  prototypical 
tissues  at  3  T  over  a  range  of  RF  phase  increments  cpQ  — 
0°-180°  (white  matter  (WM):  T1/T2  =  1000/70  ms,  diffu¬ 
sion  coefficient  D  =  0.70  x  10-3  mm2/s;  gray  matter 
(GM):  T1/T2  =  1500/100  ms,  D  =  0.80  x  10-3  mm2/s;  cer¬ 
ebrospinal  fluid  (CSF):  TJT2  =  4000/1900  ms,  D  =  3.00 
x  10-3  mm2/s).  SPGR  signals  were  simulated  using  a  = 
[3  18]°  for  three  spoiling  regimes:  weak  [TR  =  10  ms,  AG 
=  38  mT*ms/m),  strong  [TR  =  10  ms,  AG  =  110  mT*ms/ 
m),  and  nearly  complete  spoiling  [TR  =  15  ms,  AG  =  280 
mT*ms/m).  AFI  signals  were  generated  using  aAFI  =  55° 
for  weak  (TRAFI1)2  =  10/50  ms,  AG1  =  38  mT*ms/m),  in¬ 
termediate  [TRAF H  2  =  10/50  ms,  AG1  =  110  mT*ms/m), 
strong  (TRAFii,2  =  15/75  ms,  AG1  =  280  mT*ms/m),  and 
nearly  complete  spoiling  [TRAFn  2  =  30/150  ms,  AG1  = 
450  mT*ms/m).  The  data  were  then  fitted  using  the  VFA 
and  VAFI  (single  SPGR  measurement,  a  =  3°)  methods 
in  two  ways,  one  assuming  ideal  AFI  spoiling  with  noni¬ 


deal  SPGR  data,  and  one  assuming  ideal  SPGR  spoiling 
using  nonideal  AFI  data.  The  accuracy  of  T1  estimates 
was  then  assessed. 

Phantom  Experiments 

To  test  the  effects  of  the  suboptimal  AFI  regime  [TR  « 
TJ  on  both  VFA  and  VAFI,  a  set  of  10  phantoms  were 
created  from  deionized  water  doped  with  Gd-BOPTA 
to  concentrations  of  C  =  [0.00  0.10  0.25  0.50  1.00  2.00 
2.50  3.00  3.50  4.00]  mM.  The  phantoms  were  placed  in 
a  heated  water  bath  lightly  doped  with  Gd-BOPTA  and 
held  at  36.8°C  ±  0.5°C  for  the  duration  of  the  experi¬ 
ment.  SPGR  scans  were  acquired  with  TR  =  25  ms  and 
a  =  [4  10  20  30  40  50]°  and  an  AFI  scan  with  TRAFn  = 
30  ms,  aAFI  =  55°  in  the  nearly  complete  spoiling  re¬ 
gime.  Voxel  size  was  2.2  x  2.2  x  2  mm3.  All  SPGR  FAs 
were  used  in  VFA  and  VAFI  to  ensure  a  precise  mea¬ 
surement  over  the  wide  range  of  expected  T1  times. 
The  mean  and  standard  deviation  of  R1  [=1  /Ta)  were 
calculated  in  a  single  slice  region  of  interest  encom¬ 
passing  each  individual  phantom  and  used  to  fit  the 
equation: 

Ri(C)  =  Ri,o  +  ri  •  C,  [6] 

where  ra  is  the  relaxivity,  and  R1;0  is  longitudinal  relaxa¬ 
tion  rate  for  free  water.  The  linear  fit  was  weighted 
according  to  the  measured  variance  in  the  R1  values. 

It  was  shown  before  that  simultaneously  increasing 
TR  and  acceleration  factor  R  while  maintaining  a  con¬ 
stant  acquisition  time  constant  leads  to  an  improvement 
of  the  overall  SNR  efficiency  of  several  steady-state  gra¬ 
dient-echo  techniques  (38).  In  the  context  of  AFI,  in 
addition  to  signal  boost  from  longer  TR,  this  increase 
also  offers  an  opportunity  to  increase  the  length  of  the 
spoiler  gradient  and  improve  the  accuracy  of  FA  quanti¬ 
fication  (27).  To  confirm  this  effect  for  the  AFI 
sequence,  a  uniform  bottle  phantom  doped  with  Gd- 
DTPA  to  a  of  1.4  s  was  scanned  using  an  eight-chan¬ 
nel  receive-only  phased  array  knee  coil.  Two  AFI  scans 
were  acquired  (aAFI  =  55°,  nearly  complete  spoiling 
with  TR  =  30  ms,  and  weak  spoiling  with  TR  =  10  ms, 
1.5  mm3  voxel).  The  TR  =  30  ms  data  was  subsequently 
undersampled  by  a  factor  of  3  in  the  phase  encode 
direction  and  reconstructed  with  sensitivity  encoding 
(SENSE),  such  that  the  experiment  time  matched  AFI 
with  TR  =  10  ms.  FA  maps  using  weakly  spoiled  AFI 
[TR  =  10  ms)  were  also  generated  for  comparison.  Noise 
in  the  resulting  FA  maps  was  estimated  by  taking  the 
standard  deviation  in  a  small  region  of  interest  over  the 
center  of  the  bottle  to  avoid  regions  of  artifact  in  the 
poorly  spoiled  AFI  image. 

In  Vivo  Imaging 

To  demonstrate  the  application  of  VAFI  to  in  vivo  brain 
imaging,  one  healthy  human  volunteer  was  scanned. 
Informed  written  consent  was  obtained  in  accordance 
with  the  local  institutional  policy.  SPGR  scans  were 
acquired  with  TR  =  10  ms  and  aSPGR  =  [3  18]°,  which  is 
optimal  for  T1  =  1.2  s  as  described  in  (33).  To  demon¬ 
strate  the  benefits  of  increased  TR  on  spoiling,  two  AFI 
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scans  were  acquired  with  aAFI  =  55°  in  the  weak  regime 
(XT?AFI1  2  =  10/50  ms,  AG1  =38  mT*ms/m)  and  in  the 
complete  regime  (TRAFI1)2  =  30/150  ms,  AG1  =  450 
mT*ms/m).  The  last  AFI  scan  was  repeated  with  SENSE 
(37)  acceleration  factor  R  =  3  to  demonstrate  the  feasibil¬ 
ity  of  an  overall  scan  time  reduction.  All  scans  were 
acquired  with  a  128  x  96  x  44  matrix  and  an  isotropic 
resolution  of  2  x  2  x  2  mm3  placed  near  the  center  of 
the  brain.  The  scan  time  was  1:30  min  for  SPGR,  4:13 
min  for  AFI  with  TRA¥n,2  =  10/50  ms,  and  12:40  min  for 
AFI  with  TRAFI12  =  30/150  ms  and  R  =  1.  For  AFI  with 
R  =  3,  the  scan  time  was  4:13  min. 


RESULTS 

Simulations 

Figure  1  shows  the  results  of  simulations  comparing  the 
accuracy  and  precision  of  the  methods.  Over  the  range 
investigated,  the  standard  AFI  method  shows  a  system¬ 
atic  underestimation  of  FA  with  decreasing  T1  (>5%  for 
TR/T1  >  0.275).  VAFI  FA  errors  were  less  than  0.03%. 
T1  estimates  for  VFA  were  systematically  overestimated 
at  small  TR/T1}  due  to  propagation  of  AFI  errors  into 
VFA  Ta  mapping  (>5%  for  TR/T1  >  0.160).  VAFI  was 
able  to  restore  the  accuracy  of  T1  estimates  (errors  less 
than  0.07%).  TaNR  of  VAFI  exceeded  that  of  AFI  for  a 
TRIT1  <  0.085,  and  aNR  was  comparable  for  TR/T1  < 
0.150.  For  shorter  Tlt  the  AFI  method  actually  shows 
better  noise  properties  of  FA  estimates  than  VAFI;  how¬ 
ever,  this  may  be  artifactual  as  in  this  region  the  FA  and 
T1  estimates  are  also  highly  inaccurate  and  generally  not 
useful.  The  peak  TaNR  with  these  particular  settings  was 
observed  for  a  TR/T1  of  0.06;  however,  it  will  be  noted 
in  the  next  paragraph  that  7\NR  can  be  optimized  for  a 


FIG.  1.  Monte  Carlo  simulations  of  standard  VFA  (dotted  line)  and 
VAFI  (dash  line  for  two  flip  angle  fit  and  solid  line  for  one  flip  angle 
fit)  over  a  range  of  TRAF^/T^  for  flip  angle  accuracy  (a),  flip  angle- 
to-noise  ratio  (b),  T-i  accuracy  (c),  and  A-to-noise-ratio  (d). 
Results  show  a  significant  underestimation  of  FA  for  short  A 
times,  which  will  propagate  into  standard  VFA  measurements, 
causing  an  overestimation  of  Tv  VAFI  remains  accurate  over  the 
entire  range  investigated  for  both  one  and  two  SPGR  flip  angle 
options. 


cl  SPGR  FA  [degrees]  b  SPGR  FA  [degrees] 


. trafii/ti  =  0/150 - trafii/ti  =  0  030 

- TRAFI1/T1  =  0  050  - TRAFI1/T1  =  0  021 


FIG.  2.  Monte  Carlo  simulations  demonstrating  the  impact  of  flip 
angle  selection  in  the  VAFI  technique  on  the  7Vto-noise-ratio  (a, 
c)  and  flip  angle-to-noise-ratio  (b,  d)  for  a  TRAfii/Ti  =  0.150  (dot¬ 
ted  line),  0.050  (dot-dash  line),  0.030  (dash  line),  and  0.021  (solid 
line).  As  is  the  case  for  VFA,  there  exists  an  SPGR  flip  angle  which 
optimizes  precision  (a). 


specific  T1  of  interest  by  the  choice  of  an  appropriate 
SPGR  FA.  Remarkably,  VAFI  with  two  SPGR  measure¬ 
ments  shows  more  uniform  noise  performance  in  the 
range  of  TR/Tt  than  VAFI  with  single  SPGR  measure¬ 
ments,  which  is  consistent  with  improved  uniformity  of 
noise  performance  of  VFA  T1  mapping  with  extended  set 
of  SPGR  measurements  (more  than  two)  (39). 

Figure  2  shows  the  results  of  VAFI  TaNR  and  aNR  as  a 
function  of  TRIT1  and  SPGR  FA.  The  choice  of  FA  has 
an  impact  on  the  7\NR,  and  just  as  the  case  for  standard 
VFA,  there  exists  an  optimum  angle  that  maximizes  pre¬ 
cision  for  a  given  T1  time.  The  optimum  angles  for  T1  = 
200,  600,  1000,  and  1400  ms  were  7.5°,  5.5°,  4°,  and  3.5°, 
respectively  ( TR  =  10  ms).  As  long  as  the  SPGR  FA  is 
within  a  reasonable  range  about  the  optimum  value, 
there  is  no  significant  impact  on  the  precision  of  the  FA 
estimate.  The  choice  of  AFI  FA  has  considerable  impact 
on  the  precision  of  both  VAFI  T1  and  FA  measurements, 
although  in  this  case  there  is  no  single  optimum  value 
for  a  specific  T1  time;  the  precision  of  both  T1  and  FA 
increases  linearly  as  AFI  FA  increases. 

Figure  3  shows  the  results  of  simulations  assuming 
ideal  AFI  spoiling  and  nonideal  SPGR  spoiling  for  WM 
and  GM  tissues.  For  standard  VFA,  increasing  the  gradi¬ 
ent  spoiling  minimizes  the  dependence  of  T1  accuracy 
on  the  RF  phase  increment.  However,  the  choice  of  a 
proper  phase  increment  (i.e. ,  the  choice  of  cpQ  which 
crosses  or  approaches  the  zero  axis)  and  its  accuracy  and 
stability  remains  essential  for  accurate  quantification.  In 
contrast,  the  use  of  VAFI  greatly  reduces  spoiling-related 
errors  over  a  wide  range  of  phase  increments,  even  in 
the  weak  SPGR  spoiling  regime.  Specifically,  in  the 
weak  SPGR  spoiling  regime,  if  one  chooses  a  phase 
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FIG.  3.  Simulation  results  showing 
the  dependence  of  T- \  accuracy  on 
SPGR  phase  increment  for  weak  (a, 
b),  intermediate  (c,  d),  and  strong  (e, 
f)  spoiling  regimes  in  prototypical 
WM  (7V7~2  =  1000/70  ms,  D  =  0.70 
x  10-3  mm2/s)  (a,  c,  e)  and  gray 
matter  (T^/T2  =  1500/100  ms,  D 
0.80  x  10-3  mm2/s)  (b,  d,  f)  tissues 
using  standard  VFA  (dotted  line)  and 
VAFI  (dash  line  for  two  flip  angle  fit, 
solid  line  for  one  flip  angle  fit).  Simu¬ 
lations  assume  an  ideally  spoiled 
AFI  acquisition.  Results  for  CSF 
show  similar  profiles  (not  shown). 
Results  show  that  spoiling  errors 
and  dependence  on  proper  phase 
increment  are  greatly  reduced  for 
the  VAFI  technique  with  the  best 
performance  achieved  for  the  single 
flip  angle  VAFI  fit. 
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increment  that  does  not  fall  near  one  of  the  sharp  peaks 
(ex:  cp0  =  10°-35°),  VAFI  with  two  SPGR  FAs  achieves 
an  error  in  T1  estimates  of  <0.80%  (WM),  <1.10%  (GM), 
and  <1.75%  (CSF,  not  shown  in  plot)  and  an  error  in  FA 
<  0.26%  (all  tissues).  Remarkably,  these  errors  and  the 
dependency  on  the  phase  increment  are  further  reduced 
in  VAFI  with  single  SPGR  FA  (Ta  errors:  <0.54%  in 
WM,  <0.60%  in  GM,  and  <0.32%  in  CSF).  The  interme¬ 
diate  and  strong  spoiling  regimes  produce  similar 
results,  with  fewer  sharp  peaks  to  avoid. 

Figure  4  shows  the  results  of  simulations  assuming 
nonideal  AFI  spoiling  and  ideal  SPGR  spoiling  for  WM 
and  GM.  The  behavior  of  spoiling-related  errors  is  simi¬ 
lar  for  the  VFA  and  VAFI.  Unlike  the  case  of  SPGR  spoil¬ 
ing,  the  choice  of  a  proper  AFI  phase  increment  is  essen¬ 
tial  for  accurate  FA  and  T1  quantification  for  all 
techniques.  The  proper  phase  increment  (i.e.,  the  choice 
of  cp0  which  crosses  zero  axis)  depends  on  the  spoiling 
regime,  tissue  type,  and  estimation  technique  used.  For 
weak  spoiled  AFI,  the  optimal  values  of  cpQ  are  29° 
(WM),  20°  (GM),  and  50°  (CSF),  whereas  in  the  nearly 
complete  spoiling  regime,  these  values  are  39°  (WM),  41° 
(GM),  and  any  cpD  for  CSF  (not  shown  in  plot).  For  the 
VAFI  technique  in  the  weak  regime,  these  values  are  30° 
(WM),  25°  (GM),  and  55°  (CSF),  whereas  in  the  nearly 
complete  regime  these  values  are  43°  (WM),  44°  (GM), 
and  any  cpQ  for  CSF.  However,  small  deviations  in  cpQ 
about  the  optimal  value  produce  minimal  errors  in  FA 


and  Ta  accuracy  for  all  tissues  when  the  nearly  complete 
AFI  regime  is  used. 

Phantom  Experiments 

Figure  5  compares  FA  maps  in  Gd  phantoms  estimated 
using  standard  AFI  and  VAFI.  A  significant  discontinu¬ 
ity  between  the  short  T1  vials  (phantoms  6-10)  and  the 
background  media  (longer  Ta)  can  be  observed  in  the 
AFI  FA  map  (Fig.  5a).  These  errors,  which  are  due  to  a 
violation  of  the  AFI  assumption  (Eq.  1),  are  minimized 
in  the  VAFI  FA  map  (Fig.  5b),  which  agrees  well  with 
the  expected  slowly  varying  value  of  FA  in  the  imaging 
plane.  The  growing  differences  between  local  FA  values 
of  AFI  and  VAFI  with  Gd  concentration  may  be  appreci¬ 
ated  in  Table  1.  A  similar  trend  may  be  observed  in  the 
calculated  7\  values,  with  T1  errors  reaching  19%  for  the 
highest  Gd  concentration.  Corresponding  R1  measure¬ 
ments  versus  Gd  concentration  show  expected  linear 
relaxivity  curves  (Eq.  6)  for  VAFI  (Fig.  6).  The  VAFI- 
derived  relaxivity  for  the  Gd-BOPTA  solution  was  3.9361 
s-1  mM-1  [R2  =  0.9996),  which  compares  well  to  the  ref¬ 
erence  value  of  4.0  s-1  mM-1  at  37°C  (personal  commu¬ 
nication  with  Dr.  Sophie  Laurent,  NMR  and  Molecular 
Imaging  Laboratory,  University  of  Mons,  Belgium).  This 
relaxivity  measurement  was  significantly  underestimated 
by  standard  AFI  VFA  (ra  =  3.5192  s-1  mM-1,  R2  = 
0.9935). 
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FIG.  4.  Simulation  results  showing  the  dependence  of  G  and  flip  angle  accuracy  on  AFI  RF  phase  increment  for  weak  (a-d),  intermedi¬ 
ate  (e-h),  strong  (i-l),  and  nearly  complete  (m-p)  gradient  spoiling  (see  text  for  details)  using  standard  (dotted  line)  and  VAFI  (dash  line 
for  two  flip  angle  fit,  solid  line  for  one  flip  angle  fit). 


Figure  7  illustrates  the  effect  of  increasing  the  SENSE 
acceleration  factor  on  the  accuracy  and  precision  of  AFI 
maps  for  the  strong  spoiling  regime  and  a  comparison 
with  an  AFI  map  in  the  short  TR  (weak  spoiling)  regime. 
Simulated  AFI  steady-state  signals  grow  with  increasing 
TR  faster  than  the  SNR  loss  due  to  R-factor,  giving  rise 
to  an  overall  increase  in  SNR  efficiency  for  realistic 
reduction  factors  [R  <  8)  and  constant  scan  time  (Fig. 
7a).  Indeed,  long  TR  (strong  spoiling)  AFI  (Fig.  7b)  dem¬ 
onstrates  improved  precision  over  a  short  TR  (weak 
spoiling)  AFI  (Fig.  7c)  obtained  in  an  equivalent  scan 
time  [R  =  3).  The  weak  spoiling  regime  in  the  short  TR 
AFI  also  leads  to  significant  artifacts  in  the  estimated  FA 
map  (Fig.  7c). 


In  Vivo  Imaging 

Figure  8  shows  a  representative  axial  T1  map  from  the 
in  vivo  experiment  using  fully  sampled  well-spoiled 
(Fig.  8a, b)  and  poorly  spoiled  (Fig.  8d,e)  AFI  sequences. 
Table  2  shows  the  values  taken  from  region  of  interests 
positioned  in  WM,  GM,  and  CSF.  For  the  weak  AFI 
spoiling  regime  [TR1  =  10  ms),  incomplete  spoiling  leads 
to  underestimated  FA  measurements  and  overestimated 
T1  for  all  tissues,  which  is  in  accordance  with  theoretical 
predictions  (see  Fig.  4  and  Ref.  27).  In  the  complete 
spoiling  regime,  both  VFA  and  VAFI  compare  similarly 
to  previously  published  values  in  the  selected  WM/GM 


region  of  interests  (40),  while  VAFI  reaches  comparable 
or  even  better  precision  in  less  time  (using  one  less 
SPGR  measurement).  Remarkably,  in  CSF,  VAFI  pro¬ 
vided  superior  accuracy  of  Ta  estimates  based  on  the 
previously  published  reference  CSF  Ta  values  at  3  T 
(41).  The  errors  observed  for  the  standard  VFA  method 
may  be  due  to  the  fact  that  the  spoiling  regime  used  in 
this  experiment  resulted  in  incomplete  spoiling  for  the 
higher  SPGR  FA  and  nearly  complete  spoiling  for  the 
smaller  SPGR  FA.  Because  VAFI  uses  only  the  smaller 
SPGR  FA,  the  spoiling  errors  in  SPGR  are  less  likely  to 
affect  T1  estimates,  as  also  follows  from  the  theoretical 


FIG.  5.  A  map  of  flip  angle,  in  degrees,  of  the  T ^  phantom  gener¬ 
ated  using  the  standard  AFI  (a)  and  VAFI  (b)  techniques.  The  nomi¬ 
nal  FA  was  55°  (see  Table  1  for  vial  labeling).  A  significant  error  not 
related  to  the  slowly  varying  B- \  spatial  distribution  can  be  seen  in 
the  standard  case  (a),  where  a  significant  discontinuity  in  measured 
FA  exists  between  the  background  media  and  very  short  G  vials 
(#5-10).  This  error  is  corrected  using  the  VAFI  technique  (b). 
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Table  1 

Vial  Numbers,  Gd-BOPTA  Concentrations,  TA,  and  FA  Measurements  from  the  71  Phantom  (see  Fig.  7) 


Vial  # 

Gd  Cone.  (mM) 

T1  VFA  (s)  ±  Stdev 

T1  VAFI  (s) 

±  Stdev 

FA  AFI  (°) 

±  Stdev 

FA  VAFI  (°) 

±  Stdev 

1 

0.00 

4.181 

±  0.403 

4.057 

+ 

0.269 

54.1 

+ 

2.04 

54.8 

+ 

1.97 

2 

0.10 

1.523 

±  0.085 

1.494 

+ 

0.071 

61.8 

+ 

1.04 

62.0 

+ 

1.06 

3 

0.25 

0.811 

±  0.032 

0.790 

+ 

0.021 

61.6 

+ 

0.97 

62.0 

+ 

0.99 

4 

0.50 

0.451 

±  0.014 

0.432 

+ 

0.012 

49.3 

+ 

1.97 

49.7 

+ 

2.03 

5 

1.00 

0.247 

±  0.008 

0.230 

+ 

0.008 

60.1 

+ 

0.51 

61.4 

+ 

0.51 

6 

2.00 

0.137 

±  0.004 

0.122 

+ 

0.002 

60.2 

+ 

0.24 

63.1 

+ 

0.27 

7 

2.50 

0.111 

±  0.004 

0.097 

+ 

0.002 

50.9 

+ 

1.38 

53.7 

+ 

1.53 

8 

3.00 

0.093 

±  0.003 

0.083 

+ 

0.001 

48.6 

+ 

0.93 

51.8 

+ 

1.10 

9 

3.50 

0.086 

±  0.003 

0.072 

+ 

0.001 

51.9 

+ 

0.72 

56.6 

+ 

0.86 

10 

4.00 

0.079 

±  0.004 

0.064 

+ 

0.001 

46.6 

+ 

1.05 

51.2 

+ 

1.29 

analysis  (see  Fig.  3).  Finally,  we  demonstrate  a  SENSE- 
accelerated,  well-spoiled  T1  map  using  VAFI,  which 
shows  high  image  quality  (Fig.  8c)  and  T1  values  (Table 
2)  comparable  to  the  nonaccelerated  VAFI  (Fig.  8b). 

DISCUSSION 

The  VFA  SPGR  T1  mapping  continues  to  be  an  accurate 
and  efficient  method  for  mapping  T1  over  a  large  3D  vol¬ 
ume;  however,  a  large  limitation  to  this  technique  is  the 
need  for  an  accurate  estimate  of  an  actual  FA,  especially 
at  higher  field  strengths  (3  T  and  above).  Among  many 
methods  for  measuring  FA,  AFI  (25)  has  emerged  as  a 
very  attractive  choice  for  correcting  VFA  measurements 
because  of  its  high  time  efficiency  and  similar  pulse 
sequence  design.  In  spite  of  these  benefits,  AFI  often 
presents  its  own  limitations  due  to  the  need  for  a  short 
TR  to  maintain  the  validity  of  the  TRAFI12  T1  assump¬ 
tion.  The  available  tradeoffs  may  lead  to  inaccurate  FA 
and  T1  mapping  and  elongated  acquisition  times  in  some 
practical  applications.  We  presented  a  novel  technique, 
VAFI,  that  provides  a  TV-independent  FA  map  and 
simultaneously  yields  accurate  T1  values.  Briefly 
described  in  (32),  VAFI  maintains  high  accuracy  in 
regimes  that  are  problematic  for  the  standard  AFI  tech- 


FIG.  6.  Plot  of  R i  vs.  gadolinium  concentration  for  the  vial  phan¬ 
toms,  using  standard  VFA  and  VAFI.  The  standard  method  shows 
a  significant  underestimation  of  R- ,  at  short  T,  times  (due  to  errors 
in  the  AFI  FA  map  propagating  into  the  T- \  estimate)  and  the  relax- 
ivity  of  Gd-BOPTA.  [Color  figure  can  be  viewed  in  the  online  issue, 
which  is  available  at  wileyonlinelibrary.com.] 


nique  [TR  «  T V),  allowing  its  use  in  applications  involv¬ 
ing  a  long  TR  (for  proper  spoiling)  or  with  Ta-shortening 
contrast  agents  (such  as  dynamic  contrast-enhanced 
imaging).  It  was  previously  noted  that  the  VFA  tech¬ 
nique  alone  provides  no  unique  solution  for  M0,  Tlf  and 
FA  (13).  With  the  addition  of  dual-TR  AFI  data,  a  unique 
solution  can  be  found  even  with  only  a  single  SPGR  data 
point.  An  interesting  approach  to  yield  both  FA  and  T1 
values  was  proposed  in  subsequent  work,  which  general¬ 
izes  AFI  into  a  pulse  sequence  of  N  arbitrary  TR  periods 
(42).  The  practical  advantage  of  the  VAFI  approach  com¬ 
pared  to  the  multiple  TR  method  is  much  shorter  data 
acquisition  times,  which  are  minimized  when  using  a 
single  SPGR  measurement  to  augment  the  standard  dual- 
TR  AFI  scan. 

Similar  to  VAFI,  DESPOTl-HIFI  (24)  is  a  variation  on 
the  SPGR  sequence,  which  also  utilizes  simultaneous  fit¬ 
ting  to  map  FA  and  7\.  To  generate  unique  information  to 
fit  for  an  additional  FA  parameter,  DESPOTl-HIFI  utilizes 
an  inversion  pulse  prior  to  SPGR  readout.  Both  VAFI  and 
DESPOTl-HIFI  require  some  assumptions  about  RF  pulse 
properties.  In  the  case  of  DESPOTl-HIFI,  it  is  assumed 
that  the  adiabatic  inversion  pulse  is  completely  insensi¬ 
tive  to  any  parameters  affecting  the  spin  system  [B0,  Blf 
Ta),  which  requires  careful  pulse  sequence  design.  In 
VAFI,  it  is  assumed  that  the  pulse  FA  scales  linearly  with 
transmit  field.  Previous  work  (35)  has  shown  that  FA  line¬ 
arity  is  a  valid  assumption  over  the  range  of  values  herein 
investigated.  DESPOTl-HIFI  may  be  more  time  efficient 
than  VAFI  as  it  avoids  a  longer  TR.  However,  DESPOTl- 
HIFI  FA  maps  may  suffer  from  decreased  accuracy  in 
areas  with  long  T1  (e.g.,  CSF,  edema,  or  T1  “black  hole” 
lesions  in  multiple  sclerosis).  As  noted  in  (24),  these 
problems  can  be  resolved  through  the  acquisition  of  addi¬ 
tional  SPGR  FAs  and  inversion  recovery-SPGR  inversion 
times  tuned  specifically  for  the  longer  Tt,  which  would 
increase  overall  scan  time  of  DESPOTl-HIFI  accordingly. 
Another  potential  concern  is  the  reduction  of  signal  in 
DESPOTl-HIFI  due  to  magnetization  inversion.  A  detailed 
comparison  of  both  techniques,  including  their  SNR  effi¬ 
ciency,  is  a  subject  of  a  future  work. 

We  demonstrated  that,  with  properly  tuned  parameters 
(proper  choice  of  AFI  and  SPGR  FA,  phase  increment  of 
RF  spoiling,  and  spoiler  gradient  area  for  the  tissue  of  in¬ 
terest),  the  accuracy  and  precision  of  VAFI  fit  can  sur¬ 
pass  that  of  the  standard  VFA  fit  (when  used  with  AFI 
for  FA  calibration)  (Fig.  1).  This  is  true  even  when 
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FIG.  7.  a:  Effect  of  SENSE  on  the  SNR  of  the  AFI  sequence  when  simultaneously  increasing  reduction  factor  R  and  TR.  Losses  due  to 
coil  geometry  (g-factor)  were  not  considered,  b:  Flip  angle  map  acquired  with  AFI  in  strong  spoiling  regime  (long  TR)  for  SENSE  acceler¬ 
ation  factor  R  =  3.  c:  Flip  angle  map  acquired  with  AFI  in  weak  spoiling  (short  TR)  regime.  Standard  deviation  of  FA  measurements  are 
11.43  and  18.21,  respectively,  for  (b,  c),  and  scan  times  are  equivalent. 


one  SPGR  measurement  is  eliminated,  further  increasing 
the  overall  time  efficiency  of  the  technique.  This 
improvement  comes  from  the  explicit  utilization  of  least 
squares  estimation  for  solving  Eq.  9,  a  benefit  which  is 
not  available  with  linearized  calculation  of  the  FA  and 
T1  in  the  standard  AFI  and  VFA  methods.  However,  this 
improvement  is  at  the  cost  of  longer  processing  times 
than  those  required  for  original  AFI  and  VFA 
approaches.  An  alternative  and  potentially  faster  minimi¬ 
zation  approach  is  to  solve  VAFI  Eq.  5  through  inter¬ 
leaved  estimations  of  FA  by  AFI  and  T1  maps  by  VFA 
(43).  However,  this  approach  obviously  requires  a  full 
set  of  SPGR  measurements  to  allow  a  separate  VFA  fit, 
and  its  convergence  and  noise  performance  have  yet  to 
be  studied.  Image  reconstruction  can  be  easily  performed 
after  data  acquisition,  similar  to  the  majority  of  quantita¬ 
tive  parameter  mapping  methods.  The  current  implemen¬ 
tation  uses  a  brute-force  optimization  method,  and  as  a 
result  the  algorithm  requires  lengthy  processing  times  (~ 
42  milliseconds  per  voxel  or  about  6  h  for  a  high-resolu¬ 
tion  whole-brain  dataset).  An  important  next  step  will  be 
to  substantially  reduce  this  time  through  the  use  of  more 
efficient  optimization  methods  utilizing  a  C-code  imple¬ 
mentation  or  commodity  graphics  processing  unit  hard¬ 
ware  (44).  An  additional  speed  up  is  expected  from  initi¬ 


alizing  the  algorithm  using  Ta/FA  values  from  the 
regular  VFA/ AFI  fitting  procedures. 

A  major  advantage  of  the  traditional  VFA  technique  is 
the  ease  with  which  acquisition  settings  can  be  opti¬ 
mized,  either  for  a  single  T1  time  using  a  pair  of  “ideal” 
FAs  (33)  or  using  sets  of  combined  ideal  angles  to  image 
over  a  wide  range  of  T1  (10).  Our  simulation  results  show 
that  VAFI  can  be  similarly  optimized  through  the  proper 
choice  of  FA  for  the  proton-density-weighted  SPGR  mea¬ 
surement.  Because  of  the  nonlinear  nature  of  VAFI,  we 
found  this  optimal  angle  through  a  numerical  solution 
instead  of  an  empirical  formula.  It  is  important  to  note 
that  the  optimal  choice  of  this  FA  for  VAFI  is  neither  the 
Ernst  angle  nor  the  same  as  the  smallest  ideal  angle  for  a 
VFA  experiment.  However,  we  observed  that  the  optimal 
SPGR  FA  is  always  close  (within  1°)  to  the  smallest  of  the 
“ideal”  angles  optimizing  VFA  experiment  design  (33). 
Therefore,  the  “ideal”  angle  formula  may  still  be  utilized 
as  a  close  approximation  to  the  optimal  VAFI  SPGR  FA 
when  designing  experiments.  Simulations  also  showed 
that  VAFI  with  multiple  SPGR  FAs  will  result  in 
improved  efficiency  in  the  wider  T1  range,  similar  to  VFA 
T1  mapping  with  optimized  experiment  designs  (39,45). 

A  major  limitation  of  both  AFI  and  VFA  is  the  need 
for  proper  spoiling  parameters  to  ensure  accurate 


AFI  TR  =  30 
(Complete  Spoiling) 


AFI  TR  =  10 
(Weak  Spoiling) 


VFA  VAFI-1FA,  R=1  VAFI-1FA,  R=3  VFA  VAFI-1FA 


FIG.  8.  In  vivo  G  maps  in  units  of  seconds  in  a  healthy  human  volunteer  comparing  nearly  complete  AFI  spoiling  (a-c)  versus  weak 
spoiling  (d,  e)  using  standard  VFA  versus  single  flip  angle  VAFI  techniques.  For  a  well-spoiled  sequence,  both  VFA  and  VAFI  produce  G 
maps  of  comparable  quality  (see  Table  2  for  region  of  interest  measurements),  even  though  VAFI  uses  one  less  SPGR  acquisition. 
Accelerated  VAFI  demonstrates  similar  image  quality  (c).  The  result  of  a  poorly  spoiled  AFI  sequence  is  visibly  elevated  G  estimates. 
This  is  even  worse  in  the  case  of  VAFI,  where  inconsistent  signal  levels  between  the  AFI  and  SPGR  datasets  produce  visible  artifacts 
near  the  center  of  the  image  (arrow  in  e). 
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Table  2 

Region  of  Interest  Measurements  of  T ^  Comparing  Well-Spoiled  Versus  Poorly  Spoiled  AFI  Using  the  VFA  and  VAFI  Techniques, 
and  Literature  Values 


Weakly  spoiled  AFI 
(Tft-i  =  10  ms) 

Completely  spoiled  AFI  {TR^ 

=  30  ms) 

Literature 

VFA 

VAFI 

VFA 

VAFI  (R  =  1) 

VAFI  (R  =  3) 

CSF 

6.653±  1.986 

5.285  ±  1.483 

5.328  ±  1.645 

3.944  ±  0.767 

3.649  ±  0.674 

4.163  ±  0.26341 

WM  CC  genu 

1.094  ±  0.066 

1.215  ±  0.085 

0.918  ±  0.036 

0.919  ±  0.035 

0.897  ±  0.038 

0.99040 

WM  CC  splenium 

1.150  ±  0.044 

1.167  ±  0.071 

1.068  ±  0.045 

0.966  ±  0.038 

0.930  ±  0.050 

1 .02040 

WM  frontal 

1.068  ±  0.050 

1.163  ±  0.057 

0.992  ±  0.023 

1.004  ±  0.023 

0.976  ±  0.027 

0.94040 

WM  internal  cap 

1.169  ±  0.050 

1.168  ±  0.047 

1.031  ±  0.047 

0.963  ±  0.035 

0.985  ±  0.036 

0.98040 

GM  caudate 

1.818  ±  0.132 

1.982  ±  0.177 

1.505  ±  0.071 

1.535  ±  0.085 

1.480  ±  0.100 

1.41040 

GM  putamen 

1.676  ±  0.145 

1.704  ±  0.171 

1.470  ±  0.086 

1.372  ±  0.076 

1.368  ±  0.082 

1 .33040 

GM  thalamus 

1.744  ±  0.101 

1.740  ±  0.104 

1.446  ±  0.095 

1 .344  ±  0.094 

1.336  ±  0.108 

1 .22040 

quantification  (27).  Although  strong  AFI  spoiling  is  still 
a  necessity  for  VAFI,  simulations  and  in  vivo  studies 
have  shown  that  the  need  for  strong  SPGR  spoiling  in 
the  case  of  VAFI  may  be  greatly  reduced  if  a  single  FA 
SPGR  measurement  is  used.  This  may  allow  reduction  of 
spoiling  gradients  for  SPGR  part  of  VAFI  and,  as  a  result, 
more  time-efficient  overall  mapping  of  the  FA  and  T1 
than  that  is  available  with  AFI  and  VFA  alone.  Notably, 
the  dependence  of  T1  accuracy  on  the  choice  of  RF 
phase  increment  becomes  much  smaller.  This  is  similar 
to  a  previous  modification  of  the  VFA  method,  which 
proposed  an  SPGR  RF  phase  increment  near  the  plateau 
of  the  curve  rather  than  near  a  zero  crossing  in  order  to 
reduce  the  sensitivity  of  T1  to  this  value  (46).  However, 
the  previous  method  requires  an  additional  correction 
to  the  subsequent  Ta  maps,  which  is  not  necessary 
in  VAFI. 

One  must  consider  the  time  efficiency  of  VAFI  com¬ 
pared  to  the  traditional  VFA  and  AFI  approaches  when 
used  separately.  In  the  context  of  ^-insensitive  FA  map¬ 
ping,  VAFI  requires  a  short  TR  SPGR  acquisition  in  addi¬ 
tion  to  AFI,  which  may  increase  the  overall  scan  time 
insignificantly  (by  6%  in  our  unaccelerated  version  com¬ 
pared  to  standard  AFI  scan  time).  In  the  context  of  T1 
mapping,  VAFI  requires  that  both  the  AFI  and  SPGR 
scans  be  matched  in  spatial  resolution,  at  the  resolution 
of  interest  for  T1  mapping.  Although  this  eliminates  the 
need  for  one  SPGR  measurement  without  a  loss  of  preci¬ 
sion,  an  AFI  scan  with  the  optimal  n  =  TRAFi2f  TRAFn  = 
5  is  six  times  longer  than  a  normal  SPGR  acquisition, 
which  may  make  straightforward  high  resolution  (~  1 
mm3)  T1  mapping  with  VAFI  less  efficient  than  with  al¬ 
ternative  VFA-based  T1  mapping  techniques  (24). 
Although  such  full  resolution  VAFI  may  potentially 
improve  accuracy  of  FA/ TV  mapping  if  there  is  a  partial 
volume  effect  between  very  short  (problematic  AFI  re¬ 
gime)  and  long  T1  components,  several  options  may  help 
minimize  imaging  times  and  simultaneously  maintain 
accurate  FA/Ta  mapping  with  VAFI.  One  approach  is  to 
acquire  an  AFI  dataset  at  a  lower  spatial  resolution  along 
with  a  standard  VFA  acquisition  of  two  SPGR  images  at 
the  desired  resolution  and  to  use  low-resolution  VAFI  to 
fit  for  FA.  Such  method  would  still  benefit  from  high  ac¬ 
curacy  due  to  proper  spoiling  and  T1  independence  of 
FA  maps,  but  would  not  benefit  from  a  reduced  number 


of  VFA  SPGR  measurements  and  a  gain  in  T1  precision 
from  simultaneous  least  squares  fitting.  As  we  confirmed 
in  this  article,  parallel  imaging  is  also  a  very  suitable 
option  for  acceleration  of  VAFI,  as  the  signal  boost  from 
longer  AFI  TR  efficiently  compensates  SNR  losses  from 
parallel  MRI  (38),  an  effect  mostly  unavailable  for  inher¬ 
ently  short  TR  techniques  (24).  Parallel  MRI  was  used  to 
accelerate  AFI  [R  =  3)  resulting  in  the  acquisition  of  a 
VAFI  T1  mapping  in  under  5  min  with  a  minimal  effect 
on  the  map  quality  and  accuracy  (<3.6%  error  in  WM/ 
GM  Ta).  The  method  can  be  further  accelerated  for 
increased  coverage  or  higher  resolution  imaging  using 
systems  with  larger  numbers  of  receive  channels  or 
lower  resolution  VAFI  mapping  as  discussed  above.  The 
other  recommended  parameters  for  AFI  part  of  in  vivo 
protocol  are  TRAFn  =  30  ms,  TRAF I2  =  150  ms,  aAFI  = 
55°,  AG1  =  450  mT*ms/m,  cp0  =  44°,  and  a  SENSE  accel¬ 
eration  factor  of  3.  Recommended  SPGR  parameters 
are  TR  =  10  ms,  aSPGR  =  3°  AG  =  38  mT*ms/m,  and 
cp0  =  30°. 

CONCLUSION 

The  VFA  method  is  a  rapid  and  accurate  method  for 
measuring  Ta,  as  long  as  it  can  be  paired  with  an  equally 
rapid  and  accurate  method  for  measuring  and  correcting 
FA.  In  this  study,  we  presented  a  method  that  combines 
VFA  and  AFI  into  a  single  procedure,  named  VAFI.  This 
allows  the  TR  of  the  AFI  sequence  to  be  extended  to 
accommodate  a  large  spoiler  gradient  for  proper  quantifi¬ 
cation,  while  removing  the  main  AFI  requirement  of 
TRAFn,2  Ti.  We  showed  in  both  Monte  Carlo  simula¬ 
tions  and  phantom  scans  that  this  method  is  highly 
accurate  at  estimating  both  FA  and  T1  in  the  presence  of 
very  short  T1  times,  with  greater  precision  than  the  tradi¬ 
tional  VFA  approach.  We  showed  in  simulations  and  in 
vivo  scans  that  proper  spoiling  of  the  AFI  sequence  is 
essential  to  accurate  quantification,  while  the  spoiling 
requirements  of  the  SPGR  portion  of  the  experiment  are 
considerably  reduced,  resulting  in  a  much  shorter  SPGR 
acquisition  time.  The  proposed  VAFI  method  has  a  high 
potential  to  become  a  widely  used  TJ FA  mapping 
approach  in  a  variety  of  high-field  imaging  applications, 
particularly  those  with  shortened  T1  values  such  as 
dynamic  contrast-enhanced  perfusion  MRI. 
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Purpose:  Cross-relaxation  imaging  (CRI)  is  a  family  of  quanti¬ 
tative  magnetization  transfer  techniques  that  utilize  images 
obtained  with  off-resonance  saturation  and  longitudinal  relaxa¬ 
tion  rate  (fi-i)  maps  reconstructed  by  the  variable  flip  angle 
(VFA)  method.  It  was  demonstrated  recently  that  a  significant 
bias  in  an  apparent  VFA  R- \  estimation  occurs  in  macromole- 
cule-rich  tissues  due  to  magnetization  transfer  (MT)-induced 
biexponential  behavior  of  longitudinal  relaxation  of  water  pro¬ 
tons.  The  purpose  of  this  article  is  to  characterize  theoretically 
and  experimentally  the  resulting  bias  in  the  CRI  maps  and  pro¬ 
pose  methods  to  correct  it. 

Theory:  The  modified  CRI  algorithm  is  proposed,  which  cor¬ 
rects  for  such  biases  and  yields  accurate  parametric  bound 
pool  fraction  f,  cross-relaxation  rate  k,  and  R -i  maps.  Addition¬ 
ally,  an  analytical  correction  procedure  is  introduced  to  recal¬ 
culate  previously  obtained  parameter  values. 

Results:  The  systematic  errors  due  to  unaccounted  MT-induced 
biexponential  relaxation  can  be  characterized  as  an  overestimation 
of  R -i ,  f,  and  k,  with  a  relative  bias  comparable  with  the  magnitude 
of  f.  The  phantom  and  human  in  vivo  experiments  demonstrate 
that  both  proposed  modified  CRI  and  analytical  correction 
approaches  significantly  improve  the  accuracy  of  the  CRI  method. 
Conclusion:  The  accuracy  of  the  CRI  method  can  be  consider¬ 
ably  improved  by  taking  into  account  the  contribution  of  MT- 
induced  biexponential  longitudinal  relaxation  into  variable  flip 
angle  Ri  measurements.  Magn  Reson  Med  71:830-838,  2014. 
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relaxation  imaging;  variable  flip  angle  7-,  mapping 

Magnetization  transfer  (MT)  effect  is  commonly  recog¬ 
nized  as  a  source  of  important  information  about  tissue 
microstructure  due  to  its  sensitivity  to  immobile  macro- 
molecular  protons  not  detectable  by  the  conventional  MRI 
(e.g.,  those  associated  with  proteins  and  lipid  bilayers  of 
myelin,  collagen  matrix  in  cartilage,  muscle  fibers,  etc.) 
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(1,2).  Various  approaches  have  been  proposed  to  quantify 
the  MT  effect  in  vivo  including  empirical  indexes  charac¬ 
terizing  signal  saturation  due  to  MT  (1,3,4)  and  methods 
for  quantitative  mapping  of  specific  parameters  describing 
MT  within  the  two-pool  model  (2,5-15).  The  key  parame¬ 
ters  of  interest  in  this  model  are  associated  with  the  state 
of  macromolecular  protons,  which  is  characterized  by 
their  molar  fraction  (bound  pool  fraction,  f),  the  forward 
rate  constant  describing  cross-relaxation  with  water  pro¬ 
tons  [k],  and  the  transverse  relaxation  time  T%.  Particu¬ 
larly  promising  findings  were  reported  for  the  bound  pool 
fraction,  which  was  found  to  be  correlated  with  the  myelin 
content  in  neural  tissues  (16-19)  and  was  shown  to  be  ca¬ 
pable  of  tracking  age-related  changes  of  the  white  matter 
(WM)  myelination  in  animal  studies  (20). 

One  group  of  quantitative  MT  methods  termed  cross¬ 
relaxation  imaging  (CRI)  (6,10,12)  is  specifically  targeted  at 
mapping  the  parameters  /,  k,  and  Tf  or  their  subsets  in  the 
isolation  from  relaxation  properties  of  the  water  proton 
magnetization.  Collectively,  these  techniques  rely  on  the 
approximate  two-pool  pulsed  MT  formalism  (6,10)  where 
the  action  of  off-resonance  saturation  pulses  on  the  macro- 
molecular  proton  pool  is  described  by  an  effective  time-in¬ 
dependent  saturation  rate  calculated  for  a  square  pulse  with 
the  equivalent  power  and  duration.  Further,  direct  satura¬ 
tion  of  the  free  water  proton  pool  is  either  neglected  (6)  or 
estimated  within  the  stationary  approximation  sufficiently 
far  from  the  resonance  (10).  These  techniques  allow  three- 
dimensional  (3D)  acquisition  with  clinically  acceptable 
scan  times  and  resolution  based  on  a  limited  number  of 
MT -weighted  images  obtained  with  variable  off-resonance 
saturation.  The  common  feature  of  CRI  and  other  off-reso- 
nance  qMT  techniques  is  the  need  for  complementary  TV 
mapping,  which  allows  decoupling  of  the  two-pool  model 
parameters  from  the  longitudinal  relaxation  rate  7?a  =  l/7\ 
(5-7,9,10,12).  Because  of  the  need  for  the  fast  3D  acquisi¬ 
tion,  T1  maps  are  often  generated  using  the  variable  flip 
angle  (VFA)  method  with  a  spoiled  gradient-echo  (SPGR) 
sequence  (21).  It  has  been  realized  recently  that  the  VFA 
method  originally  derived  for  the  single-pool  model  inac¬ 
curately  describes  the  SPGR  signal,  which  often  leads  to  bi¬ 
ased  7?!  estimation  in  tissues  with  rich  macromolecular 
content.  This  bias  is  due  to  unaccounted  cross-relaxation 
caused  by  MT  effect  between  macromolecular  and  water 
protons,  which  leads  to  the  biexponential  behavior  of  longi¬ 
tudinal  relaxation  of  water  proton  magnetization  in  such  tis¬ 
sues  (22).  The  cross-relaxation  contribution  was  shown  to 
dominate  7?a  of  water  in  hydrated  collagen  (22)  and  more 
recently  in  neural  tissues  (23).  Deviation  of  longitudinal 
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relaxation  from  a  single-exponential  behavior  due  to  cross¬ 
relaxation  significantly  affects  signal  intensities  a  in  the  fast 
gradient-echo  sequences  (24,25)  and  introduces  a  bias  in 
VFA  f?!  measurements  that  depends  on  the  pulse  sequence 
parameters  and  may  reach  up  to  14-15%  in  brain  WM  (24). 

Apart  from  biasing  VFA  R1  estimation,  systematic  errors 
from  unaccounted  biexponential  relaxation  may  further 
propagate  into  the  two-pool  model  parameters  estimated  by 
quantitative  MT  techniques  such  as  CRI.  Conversely,  an 
unbiased  determination  of  R1  values  in  tissues  requires  the 
knowledge  about  the  rest  of  parameters  of  the  two-pool 
model,  which  may  not  be  practical  for  data  correction  within 
the  VFA  method  alone  (24).  In  this  study,  we  propose  the 
unified  treatment  of  VFA  and  MT  SPGR  signals  using  a 
modified  CRI  (mCRI)  analysis,  which  enables  simultaneous 
corrections  of  R1  and  two-pool  MT  model  parameters.  In  par¬ 
ticular,  we  theoretically  and  experimentally  characterize 
systematic  errors  in  CRI  caused  by  biexponential  relaxation 
in  VFA  R1  mapping  and  demonstrate  a  new  processing  algo¬ 
rithm,  which  corrects  for  such  errors  and  yields  accurate 
parametric  /,  k,  T%,  and  R1  maps.  Additionally,  we  propose 
an  analytical  correction  procedure  allowing  recalculation  of 
previously  obtained  cross-relaxation  parameter  values  with 
acceptable  residual  errors. 

THEORY 

Analytical  Theory  of  the  SPGR  Signal  in  the  Presence  of 
Cross-Relaxation 

To  estimate  the  effect  of  cross-relaxation  on  the  apparent 
R1  measured  by  the  VFA  method,  we  employ  the  previ¬ 
ously  described  pulsed  steady-state  formalism  (6,10)  in 
the  simplified  form,  where  no  off-resonance  saturation  is 
applied.  The  matrix  model  of  the  longitudinal  magnet¬ 
ization  (based  on  Eq.  [1]  in  Ref.  10)  can  be  rewritten  as: 

Mz  =  (I  — EC)_1(I— E)Me(?,  [1] 

where  Mz=  [Mf  Mf] T  corresponds  to  the  longitudinal  mag¬ 
netization  before  the  excitation  pulse,  Mf ,  Mf  are  the  longi¬ 
tudinal  magnetizations  of  free  and  bound  protons, 
respectively,  Meq  =M0[1—  /  f]T  is  the  vector  of  equilibrium 
magnetization,  I  is  the  identity  matrix,  the  matrix  E= 
exp  (RTr)  describes  relaxation  during  repetition  time  TR, 
and  the  diagonal  matrix  C=diag  (Sf,  S*,),  Sf =cos  a,  S^  =  l 
corresponds  to  the  instant  rotation  of  the  magnetization  Mf 
by  an  excitation  pulse  with  the  flip  angle  a.  The  relaxation 
matrix  R  is  defined  as  follows: 


k  -RB-kML 


[2] 


where  R ^  and  Rf  are  the  longitudinal  relaxation  rates  of 
the  free  and  bound  pool,  respectively.  Applying  the  first- 
order  approximation  to  the  exponential  terms,  Eq.  [1] 
can  be  simplified  to 


Mz  R3  (RTfl  +  In  C)  1RMeq  Tr  [3] 


S 


M0  (1  -/)  (RB  RB +RB  Jc  (1  -/)  //+ RB  k)  sin  aexp  ( -  TE/T * ) 
R^+R^k(l-f)/f+RB+k-{RB+k(l-f)/f)T^ln(cosix) 


[4] 


We  further  assume  that  Rb=Rb=Ri  similar  to  earlier 
studies  (6,12)  that  corresponds  to  the  fast  exchange  con¬ 
ditions  (26).  With  these  assumptions,  the  signal  intensity 
is  expressed  as: 

M0(l-/)f?1TflSmaexp  (-Tjj/T*) 

\  |  Mn (cos a) 

RiTr  In  (cos  a)  +  Rl+k/f 


Finally,  the  relaxation  rate  R1  can  be  neglected  in  the 
sum  with  a  much  larger  term  k/f,  thus  providing  the  sig¬ 
nal  equation 

S  *  Mo<1~/)H,r,-(fmin(cos»)Slll"Xp  (-7i/T;) 

[6] 


By  comparing  Eq.  [6]  with  the  first-order  approxima¬ 
tion  of  the  Ernst  equation  (27)  (not  shown  for  brevity), 
the  relationship  between  the  true  R1  and  its  apparent 
value  R*pp  estimated  from  VFA  data  becomes  evident: 

R^^R1/(l-f)  [7] 

Equation  [7]  demonstrates  that  the  bias  in  R1  caused 
by  cross-relaxation  is  on  the  order  of  /  and  independent 
of  the  sequence  parameters  ( TR  and  flip  angles)  if  the 
first-order  approximation  is  justified  by  a  short  TR.  More 
general  treatment  for  arbitrary  sequence  parameters  can 
be  found  in  Ref.  24  and  25. 


Effect  of  Apparent  on  the  CRI  Parameters 

The  next  goal  is  to  estimate  systematic  errors  in  the  two- 
pool  model  parameters  due  to  R^pp  in  the  CRI  method. 
We  start  our  analysis  with  a  simplified  CRI  model  (6), 
which  allows  analytical  investigation  of  errors  in  the 
estimated  parameters  in  the  regime  where  the  direct  sat¬ 
uration  effect  is  negligible  (A  >2.5  kHz)  (28).  According 
to  this  model,  the  ratio  of  signals  with  and  without  off- 
resonance  saturation  in  a  pulsed  MT  experiment  can  be 
expressed  as 


M?(Wb)  _P+(Q-1)sWb 
MB{WB= 0)  ~  P+QsWB 

where 

P=f~ 1  (fli ■ -(!-/) TRHn  (cos  a))  [9] 

Q=k_1(R1-TR1  In  (cos  a)) +  1  [10] 


Corresponding  approximated  expression  for  the 
observed  signal  can  be  explicitly  written  as  follows: 


and  WB  is  the  saturation  rate  of  the  bound  pool  scaled 
by  the  duty  cycle  of  the  saturation  pulse  [s  =  tmt/TR, 
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where  tmt  is  the  pulse  duration).  The  saturation  rate  W6 
depends  on  the  offset  frequency  and  the  flip  angle  of  the 
saturation  pulse  and  is  determined  by  the  parameter  T2 
based  on  an  appropriate  spectral  line  shape  model  [typi¬ 
cally  SuperLorentzian  (29)].  If  the  data  are  fitted  in  the 
form  given  by  Eq.  [8],  the  fitting  algorithm  searches  for 
the  optimal  parameters  P,  Q,  and  T2  regardless  of  an 
actual  i?a. 

The  effect  of  R1  errors  on  the  two-pool  model  parameters 
originates  from  the  dependence  of  the  coefficients  P  and  Q 
on  R1  (Eqs.  [9]  and  [10]).  Accordingly,  if  R*pp  is  supplied 
as  input  information,  apparent  values  /app  and  will  be 
obtained.  The  relationships  between  the  true  /  and  k  and 
their  apparent  values  can  be  derived  from  the  conditions 

P^app  ?  y aPP)  =  p[R^  J)  and  Q(R apP  ?  j^PP)  =  ^  which 

result  in  the  simple  recalculation  formulas: 

f=Cf\  k=Ck*vv  [11] 

where  the  correction  coefficient  C  is  the  same  for  both  / 
and  k  and  depends  on  TR  and  the  flip  angle: 

C=  (R*pp  -Tn1  In  (cos  a))  /  (R*pp  -TRHn  (cos  a)  +/appl?app) 

[12] 

The  above  analysis  also  suggests  that  R1  errors  should 
not  affect  T2 ,  since  W5  is  independent  of  R1  (6).  We  refer 
below  to  the  recalculation  formulas  (Eqs.  [11]  and  [12]) 
as  CRI  with  the  first-order  correction. 

Standard  CRI  and  mCRI  Approaches 

In  the  original  CRI  approach  (10),  reconstruction  of  para¬ 
metric  maps  is  performed  in  two  stages  (Fig.  1).  During 
the  first  stage,  R1  map  is  calculated  by  fitting  the  Ernst 
equation  to  VFA  SPGR  data  (21).  During  the  second  step, 
the  matrix  model  of  the  pulsed  MT  is  fitted  to  MT- 
weighted  data,  while  R1  is  supplied  as  an  external  param¬ 
eter.  Prior  to  the  fit,  the  MT  data  are  normalized  pixel- 


wise  by  an  arbitrary  reference  image  typically  obtained 
without  saturation.  The  purpose  of  this  normalization  is 
to  reduce  the  number  of  free  parameters  in  the  fit  by 
excluding  a  multiplicative  factor  commonly  referred  to  as 
proton  density,  which  absorbs  the  effects  of  spin  concen¬ 
tration,  T2*  decay,  and  the  coil  sensitivity.  The  normal¬ 
ized  data  are  then  fitted  to  the  corresponding  ratio  of 
analytical  signal  expressions  with  and  without  saturation 
Mf  (VES)/Mf  (VP5  =  o).  As  an  alternative  to  the  standard 
CRI  approach,  we  propose  the  mCRI  reconstruction  algo¬ 
rithm  to  perform  a  global  fit  of  VFA  and  MT  SPGR  data 
simultaneously  (Fig.  1).  On  the  first  stage,  similarly  to  the 
standard  CRI  approach,  the  normalization  procedure  is 
uniformly  applied  to  both  VFA  and  MT  data  to  exclude 
the  common  proton  density  term.  Then,  the  normalized 
images  are  simultaneously  fitted  to  the  corresponding 
ratios  of  analytical  signal  expressions  with  and  without 
saturation  using  Eq.  [1]  from  this  paper  and  Eq.  [1]  from 
Ref.  10,  respectively.  Accordingly,  the  global  fit  simulta¬ 
neously  yields  all  the  remaining  parameters  [Rlf  /,  k,  T%) 
and  automatically  takes  into  account  the  cross-relaxation 
contribution  into  VFA  data  (Fig.  1). 

METHODS 

Simulations 

To  estimate  errors  in  the  two-pool  model  parameters 
caused  by  the  unaccounted  biexponential  relaxation, 
synthetic  VFA  and  MT  SPGR  signal  intensities  were  gen¬ 
erated  using  Eq.  [1]  from  this  paper  and  Eq.  [1]  from  Ref. 
10,  respectively.  Then,  datasets  were  fitted  using 
the  original  CRI  algorithm.  The  first-order  correction  for¬ 
mulas  (Eqs.  [7],  [11],  and  [12])  were  also  evaluated.  Spe¬ 
cific  parameters  used  in  the  simulations  were  taken  from 
regions-of-interest  (ROI)  measurements  in  the  genu  of 
corpus  callosum  (WM)  and  thalamus  (GM)  from  full 
mCRI  fit  of  the  in  vivo  data  presented  below:  1)  WM: 
R1  =  0.97  s-1;  /=  15.36  %;  k=  2.70s-1;  Tf  =  9.84  ps;  2) 
GM:  R1  =  0.71  s-1;  /=  8.8%;  k=  1.57  s_1;  Tf  =  10.21  ps. 
To  study  the  effect  of  approximations  of  the  first  order 
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FIG.  1.  Original  CRI  and  mCRI  processing  pipelines. 
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correction  in  the  realistic  imaging  regime,  pulse 
sequence  parameters  and  sampling  scheme  in  these  sim¬ 
ulations  corresponded  to  those  used  in  the  in  vivo 
experiments  detailed  below. 

Phantom  Preparation 

Five  crosslinked  BSA  samples  (98%  bovine  serum  albu¬ 
min,  essentially  fatty  acid  free,  Sigma-Aldrich  Corp.,  St. 
Louis,  MO)  were  prepared  with  BSA  percent  weight  of 
10,  15,  20,  25,  and  30  as  described  in  Ref.  30.  The  BSA 
was  dissolved  in  distilled  water  and  was  placed  in  the 
ice  bath  for  10  min.  Then,  50  pL/mL  of  an  ice-cold  25% 
glutaraldehyde  solution  (Sigma-Aldrich  Corp.,  St.  Louis, 
MO)  were  added  to  the  BSA  solution  while  stirring  with 
a  syringe  needle.  The  samples  were  kept  at  the  room 
temperature  for  two  hours  and  then  stored  at  4°C. 

Data  Acquisition 

Imaging  experiments  were  carried  out  on  a  3.0  T  GE  Dis¬ 
covery  MR750  (GE  Healthcare;  Waukesha,  WI)  using  ei¬ 
ther  an  eight-channel  transmit/receive  knee  coil  (for 
phantom  scans)  or  eight-channel  phased  array  head  coil 
(for  a  volunteer  scan).  All  data  were  acquired  with  the 
3D  MT-weighted  SPGR  sequence  in  a  strong  spoiling  re¬ 
gime  (31)  (spoiling  gradient  area  AG  =  450  mT-ms/m,  RF 
phase  increment  169°).  In  the  phantom  experiment,  eight 
Z-spectroscopic  datasets  were  acquired  [TR/TE  =  37/2.3 
ms,  excitation  flip  angle  a  =15°)  with  the  18  ms  Fermi 
saturation  pulse  applied  at  the  offset  frequencies  A  =  2.5, 
10,  18,  26  kHz  with  two  nominal  saturation  flip  angles 
aMT  =  850°  and  1300°.  Additionally,  four  VFA  SPGR  data¬ 
sets  were  acquired  using  the  same  sequence  with  A  =  250 
kHz  to  ensure  that  the  transmitter  operates  with  identical 
gain  settings  (no  MT  effect  is  observed  at  this  frequency) 
and  flip  angles  optimized  for  the  range  of  T1  values  in  the 
phantoms  (a  =  6°,  15°,  35°,  and  50°).  The  SPGR  dataset 
with  the  highest  signal -to-noise  ratio  (a  =  15°)  was  used  as 
a  reference  image  to  normalize  both  Z-spectroscopic  and 
VFA  data  as  explained  in  the  previous  section.  All  data¬ 
sets  were  acquired  with  field  of  view  =140  X  105  X  48 
mm3  and  matrix  =  128  X  96  X  24.  Single-slice  2D  inver¬ 
sion-prepared  spin-echo  (IR)  data  were  collected  to  deter¬ 
mine  reference  T1  values  in  the  phantoms  ( TR/TE=  5000/ 
8.2  ms,  Ti  =  0.3,  0.5,  0.7,  0.9,  1.2,  1.6,  and  2  s). 

Informed  written  consent  was  obtained  from  a  healthy 
volunteer  in  accordance  with  the  local  institutional  pol¬ 
icy.  Eight  Z-spectroscopic  datasets  were  acquired  using 
the  same  pulse  sequences  (A  =  2.5,  5,  9,  13  kHz, 
o(mt  =  500°,  1100°,  a  =  10°).  Four  VFA  datasets  were 
acquired  with  flip  angles  a  =  5°,  10°,  20°,  and  30°.  The 
image  with  a  =  10°  was  used  as  a  reference  dataset.  All 
Z-spectroscopic  and  VFA  data  were  acquired  with  TR/ 
Te  =  40/2.0  ms,  field  of  view  =  240  X  180  X  80  mm3  and 
matrix  =  128  X  96  X  42.  The  total  scan  time  for  this  pro¬ 
tocol  was  35  min. 

B0  and  B1  maps  were  used  to  correct  flip  angle  and 
local  off-resonance  frequency  values  in  both  phantom 
and  volunteer  studies.  B1  maps  were  acquired  using  the 
actual  flip  angle  imaging  method  (32)  (TRtl/TRt2/TE  =  37/ 
185/2.3  ms,  a  =  55°,  field  of  view  =  240  X  180  X  80  mm, 
matrix  =  96  X  72  X  28)  with  the  strong  spoiling  regime 


(31).  B0  maps  were  calculated  from  the  3D  fat-water  sep¬ 
aration  method  known  as  “iterative  decomposition  of 
water  and  fat  with  echo  asymmetry  and  least  squares 
estimation”  with  field  of  view  =240  X  180  X  80  mm 
and  matrix  =  256  X  256  X  42  (33).  The  scan  time  for 
these  additional  B0  and  B1  mapping  sequences  was  10 
min. 

Image  Processing  and  Analysis 

The  standard  CRI  processing  workflow  (Fig.  1)  was  imple¬ 
mented  according  to  (10,34).  In  mCRI,  all  parametric  maps 
were  generated  by  fitting  the  normalized  VFA  and  MT  data 
simultaneously  as  described  above  (Fig.  1)  using  in-house- 
written  C  and  MATLAB  (MathWorks,  Natick,  MA)  software 
utilizing  a  standard  “ lsqnonlin”  function  for  nonlinear 
least  squares  voxel-based  fitting  (http://www. 

medphysics.wisc.edu/~samsonov/qmap).  All  quantitative 
parameters  estimated  in  the  phantom  scans  and  in  vivo 
were  obtained  from  parametric  maps  in  the  manually  drawn 
ROIs. 

To  study  the  association  between  BSA  content  and  the 
bound  pool  fraction  determined  by  different  processing 
approaches,  we  performed  linear  regression  analysis 
using  the  following  model: 

f=P1  *£SA%+/0  [13] 

The  strength  of  anticipated  linear  relationships  was 

assessed  using  the  Pearson’s  correlation  coefficient.  We 

additionally  tested  if  the  intercept  /0  is  significantly  dif¬ 
ferent  from  zero  (the  anticipated  value  for  /  at  BSA%  =  0). 
Statistical  significance  of  differences  between  parameter 
measurements  by  CRI,  mCRI,  and  CRI  with  the  first-order 
correction  was  assessed  using  a  paired  two-tailed  f-test. 

RESULTS 

Simulations 

Table  1  shows  the  systematic  errors  in  the  estimated  pa¬ 
rameters  in  WM  and  GM  obtained  from  the  original  CRI 
method,  mCRI,  and  CRI  with  the  first-order  correction.  As 
predicted  by  the  theory,  the  R1  bias  propagates  into  / and  k 
measurements  but  does  not  affect  Tf ,  and  it  is  most  promi¬ 
nent  for  WM.  A  major  portion  of  the  biases  are  removed  by 
the  proposed  first-order  correction,  with  R1  and  /  being 
corrected  most  efficiently  (<1%  residual  error). 

Figure  2  further  illustrates  the  dependence  of  the 
biases  on  the  fraction  of  bound  protons  in  WM.  Again, 
the  unaccounted  biexponential  relaxation  in  the  SPGR 
signal  model  introduces  a  substantial  bias  into  all  quan¬ 
titative  MT  parameters  except  for  T%  when  estimated  by 

Table  1 

Simulated  Relative  Bias  in  Two-Pool  MT  Model  Parameters 
Obtained  by  the  Original  CRI  Fit  and  Its  First-Order  Correction  due 
to  the  Unaccounted  Biexponential  Relaxation  for  WM  and  GM 


Original  CRI 

First-order  correction 

Hi  (%) 

f(%) 

k  (%> 

7?  (%) 

fli  (%) 

f(%) 

k  (%) 

WM 

16.92 

12.42 

15.20 

-0.08 

-0.96 

-0.42 

2.04 

GM 

8.88 

5.95 

6.46 

-0.07 

-0.64 

-0.24 

0.24 
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FIG.  2.  Relative  bias  in  parameters  estimated  using  original  CRI 
before  and  after  first-order  correction  vs.  bound  pool  fraction.  The 
other  two-pool  model  parameters  were  fixed  with  values  corre¬ 
sponding  to  WM. 

the  original  CRI  method.  These  errors  affect  most  signifi¬ 
cantly  ib  values,  then  k  and  /  values.  They  are  steadily 
growing  with  the  bound  pool  fraction  as  predicted  by 
Eqs.  [7]  and  [11].  The  first-order  correction  (Eqs.  [11]  and 
[12])  removes  most  of  the  bias  for  a  wide  range  of  /val¬ 
ues,  while  the  level  of  residual  errors  increases  with  an 
increase  in  /.  The  first-order  correction  is  less  accurate 
for  k  values  where  the  residual  bias  may  reach  up  to 
~5%  for  tissues  with  a  higher  macromolecular  content. 


Phantom  Studies 

Figure  3  shows  estimates  of  the  two-pool  MT  model  pa¬ 
rameters  in  the  BSA  phantoms  obtained  by  the  original 
CRI,  mCRI,  and  CRI  with  the  first-order  correction,  as 
well  as  fb  values  from  the  reference  IR  experiment.  Both 
mCRI  fit  and  the  first-order  correction  of  CRI  fit  provide 
a  good  agreement  with  the  reference  IR  fb  measure¬ 
ments,  while  the  standard  VFA  method  results  in  sub¬ 
stantial  discrepancy  with  IR  fb  estimates.  The  difference 
between  fb  values  measured  by  the  VFA  method  and  IR 
method  increases  with  BSA  concentration  as  predicted 
by  Eq.  [7]  and  the  simulations  (Fig.  2). 

Similarly,  the  original  CRI  tends  to  overestimate  bound 
pool  fraction  /  and  cross-relaxation  rate  k  in  the  phan¬ 
toms  with  a  higher  concentration  of  BSA.  Although  the 
first-order  correction  of /provides  a  good  agreement  with 
the  mCRI  fit  for  all  concentrations,  an  increase  in  the 
macromolecular  content  leads  to  a  poorer  first-order  cor¬ 
rection  for  k,  as  predicted  by  the  simulations  (Fig.  2).  T. f 
is  consistent  among  different  concentrations  of  BSA. 

Figure  4  and  Table  2  show  results  of  a  linear  regres¬ 
sion  analysis  of  /  from  phantom  data.  The  bound  pool 
fraction  strongly  correlates  with  BSA  concentration  for 
both  CRI  and  mCRI,  though  the  original  CRI  estimation 
lead  to  a  somewhat  weaker  correlation  (Table  2).  The 
intercept  of  the  fitted  line  for  mCRI  and  the  first-order 
correction  is  not  significantly  different  from  zero,  thus 
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FIG.  3.  Comparison  of  qMT  parameters  estimated  using  original  CRI  (without  and  with  first-order  correction)  and  mCRI  in  phantoms 
with  different  BSA  concentrations,  a:  R -i  measurements.  IR  R -i  values  are  shown  for  reference,  b-d:  f,  k,  and  Tf  measurements,  respec¬ 
tively.  Error  bars  indicate  standard  deviations  in  ROI  measurements  taken  from  parameter  maps. 
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FIG.  4.  Results  of  linear  regression  of  bound  pool  fraction  values 
derived  by  CRI,  CRI  with  the  first-order  correction,  and  mCRI  vs. 
BSA  concentration.  Error  bars  indicate  standard  deviations  in  ROI 
measurements  taken  from  parameter  maps. 


reflecting  the  anticipated  absence  of  the  bound  proton 
pool  at  0%  BSA  concentration.  At  the  same  time,  the 
intercept  of  the  fitted  line  for  original  CRI  is  significantly 
different  from  zero  and  appears  in  the  physically  unreal¬ 
istic  range  (~— 2.5%). 

In  Vivo  Results 

Figure  5  compares  parametric  maps  obtained  using  the 
mCRI  processing  approach  with  original  CRI  and  its 
first-order  correction.  Table  3  provides  quantitative  com¬ 
parisons  of  qMT  parameters  in  several  WM  and  GM 
regions  of  interest.  Processing  with  mCRI  resulted  in 
structurally  similar,  but  quantitatively  different  paramet¬ 
ric  maps.  For  all  parameters  except  for  T%,  the  error 
images  reveal  dependence  of  the  bias  on  the  macromo- 
lecular  content  as  predicted  by  simulations  and  phantom 
studies.  First-order  correction  efficiently  minimizes  the 
errors  in  /  and  R1.  However,  there  is  still  a  visible  resid¬ 
ual  error  for  k  values  after  the  first-order  correction,  with 
most  appreciable  residue  in  WM  regions,  which  is  in 
agreement  with  simulations  and  phantom  experiments 
(Figs.  2  and  3,  Table  1).  Statistical  analysis  of  ROI  data 
detected  significant  differences  between  CRI  and  mCRI 
for  all  parameters  except  for  T%  in  gray  matter  (Table  4). 
The  effect  of  the  first-order  correction  was  also  signifi¬ 
cant  for  R1,  f,  and  k.  There  were  significant  differences 
between  mCRI  and  the  first-order  correction  for  R1  and 
k.  These  systematic  errors  as  well  as  bias  between  mCRI 
and  the  first-order  correction  (Table  4)  are  consistent 
with  the  estimates  predicted  in  simulations  (Table  1). 


a  b 


FIG.  5.  Parametric  MT  maps  estimated  using  original  CRI  and 
mCRI  (a),  and  errors  of  the  original  CRI  and  its  first-order  correction 
[CRI  (1st)]  presented  as  difference  maps  with  respect  to  mCRI  (b). 
Note  the  consistency  between  levels  of  error  in  individual  maps  vs. 
macromolecular  content  as  revealed  by  the  bound  pool  fraction. 

DISCUSSION 

MT  imaging  offers  unique  sensitivity  to  macromolecular 
content  in  tissues  that  cannot  be  assessed  with  conven¬ 
tional  MRI.  Accurate  modeling  of  MT  processes  using 
the  two-pool  model  is  essential  to  yield  specific  parame¬ 
ters  characterizing  the  macromolecular  fraction  in  both 
normal  and  pathological  conditions.  This  study  demon¬ 
strates  on  the  example  of  CRI,  one  of  qMT  methods 
based  on  SPGR  pulse  sequence,  that  separate  treatment 
of  VFA  and  MT  data  may  cause  non-negligible  system¬ 
atic  errors  in  both  R1  and  quantitative  MT  parameters 
such  as  bound  pool  fraction  /  and  cross-relaxation  rate  k. 
Specifically,  this  type  of  errors  can  be  characterized  as 
an  overestimation  of  Rlf  f,  and  k  with  a  relative  bias 
comparable  with  magnitude  of /(Fig.  2,  Tables  1,  3,  and 
4).  As  the  effect  of  the  MT-induced  biexponential  longi¬ 
tudinal  relaxation  on  the  parameters  estimation  strongly 
depends  on  the  macromolecular  content  in  tissues  (Fig. 
2),  the  original  CRI  method  and  other  qMT  analysis  tech¬ 
niques  that  do  not  take  into  account  this  effect  may 
make  interpretation  of  quantitative  parameters 


Table  2 

Results  of  Linear  Regression  Analysis  for  Bound  Pool  Fraction  f  in  BSA  Phantoms 


Po  (P- value) 

ft 

R2 

CRI 

-2.46e-2  ±  0.18e-3  (<1(T4)* 

1.02e-2±8.63e-5 

0.965 

First-order  correction 

-0.25e-2±1.31e-3  (0.053) 

0.81e-2±6.16e-5 

0.972 

mCRI 

-0.14e-2±1.27e-3  (0.256) 

0.80e-2  ±  6.00e-5 

0.972 

Asterisk  (*)  indicated  that  the  difference  between  the  intercept  of  the  fitted  line  and  zero  is  statistically  significant. 


836 


Mossahebi  et  al. 


Table  3 

In  Vivo  Measurements  of  the  Two-Pool  MT  Model  Parameters  Using  mCRI,  Original  CRI,  and  CRI  with  First-Order  Correction 


mCRI 

Original  CRI 

First-order  correction 

Regions-of-interest 

Ri  (s'1) 

f(%) 

k  (s-1) 

7f  (ns) 

(s-1) 

f(%) 

k  (s-1) 

7f  (ns) 

«1  (s-1) 

f(%) 

k  (s-1) 

White  matter 

Cerebral  peduncle 

0.933 

14.81 

2.064 

11.560 

1.088 

16.54 

2.353 

11.556 

0.928 

14.73 

2.097 

Corona  radiata 

0.894 

13.87 

1.979 

11.731 

1.032 

15.44 

2.247 

11.722 

0.889 

13.87 

2.019 

Corpus  callosum,  genu 

0.966 

15.36 

2.702 

9.844 

1.133 

17.28 

3.115 

9.836 

0.959 

15.30 

2.758 

Corpus  callosum,  splenium 

0.942 

14.92 

2.517 

10.197 

1.105 

16.75 

2.890 

10.187 

0.940 

14.90 

2.570 

Frontal  white  matter 

0.898 

13.93 

2.161 

11.231 

1.036 

15.49 

2.445 

11.223 

0.892 

13.90 

2.196 

Internal  capsule,  anterior  limb 

0.932 

13.88 

2.327 

10.171 

1.085 

15.52 

2.651 

10.162 

0.934 

13.91 

2.378 

Internal  capsule,  posterior  limb 

0.903 

15.33 

2.048 

11.836 

1.057 

17.20 

2.356 

11.829 

0.896 

15.26 

2.092 

Middle  cerebrallar  peduncle 

0.855 

14.28 

1.993 

10.473 

0.999 

15.99 

2.220 

10.443 

0.855 

14.33 

1.990 

Occipital  white  matter 

0.924 

14.79 

1.898 

10.628 

1.080 

16.73 

2.198 

10.621 

0.919 

14.89 

1.956 

Gray  matter 

Cerebral  cortex 

0.739 

6.6 

1.102 

9.023 

0.781 

6.89 

1.122 

8.987 

0.729 

6.58 

1.072 

Caudate  nucleus 

0.655 

6.42 

1.321 

9.628 

0.702 

6.71 

1.380 

9.619 

0.657 

6.43 

1.322 

Putamen 

0.710 

7.29 

1.432 

9.736 

0.770 

7.66 

1.512 

9.728 

0.714 

7.29 

1.439 

Substantia  nigra 

0.788 

8.59 

1.371 

9.817 

0.863 

9.14 

1.462 

9.808 

0.789 

8.59 

1.375 

Thalamus 

0.707 

8.77 

1.569 

10.208 

0.780 

9.36 

1.682 

10.198 

0.711 

8.8 

1.583 

inconsistent  for  different  tissues  types  in  different  patho¬ 
logical  conditions.  It  is  also  important  to  emphasize  that 
apparent  R1  values  obtained  from  the  standard  single¬ 
compartment  VFA  technique  are  biased  proportionally  to 
/,  and  hence  represent  nonspecific  measures  containing 
contributions  from  both  the  intrinsic  T1  relaxation  prop¬ 
erties  of  water  protons  in  tissues  and  the  cross-relaxation 
between  two  proton  pools. 

We  demonstrated  that  the  accuracy  of  CRI  estimates  can 
be  noticeably  improved  with  the  proposed  combined  data 
fit  algorithm.  Our  mCRI  data  processing  approach  does 
not  require  any  additional  measurements,  thus  maintain¬ 
ing  the  same  time-efficiency  as  the  original  CRI  technique. 
The  mCRI  method  yields  unbiased  VFA  R1  estimation, 
resulting  in  a  decrease  in  R1  values  by  approximately  15% 
in  WM  and  8%  in  GM  (Tables  3  and  4),  which  agrees  well 
with  the  bias  predicted  earlier  (24).  In  addition  to  the 
mCRI  fit,  we  have  proposed  a  simple  analytical  first-order 
approximated  correction  formulas  allowing  recalculation 
of  original  CRI-based  parameters.  The  residual  bias  of  such 
correction  was  shown  both  theoretically  and  experimen¬ 
tally  to  increase  with  macromolecular  content  due  to  viola¬ 
tion  of  the  first-order  approximation.  It  may  reach  about 
2%  for  k  and  1%  for  both /and  R1  in  the  ranges  of  physio¬ 
logically  reasonable  values  for  brain  tissues.  Nevertheless, 
the  first-order  correction  employed  as  a  simple  postpro¬ 
cessing  step  in  the  original  CRI  processing  pipeline  (Fig.  1) 
may  be  a  valid  alternative  to  mCRI  when  either  faster  proc¬ 
essing  is  required  or  previously  computed  parameter  values 

Table  4 

Average  Relative  Bias  Between  mCRI  and  CRI  Without  and  With 
First-Order  Correction  Calculated  from  In  Vivo  ROI  Measurements 
(Table  3) 


CRI 

First-order  correction 

(%) 

f  (%)  k  (%) 

7f  (%) 

Ri  (%> 

f(%) 

k  (%> 

WM 

15.33 

11.34  13.21 

0.09 

-0.41 

-0.06 

1.85 

GM 

7.94 

5.39  5.2 

0.15 

0.03 

0.05 

0.06 

Bold  font  indicates  significantly  different  results  (P<0.05). 


need  to  be  refined.  As  the  residual  error  after  first-order  cor¬ 
rection  increases  with  /  care  must  be  taken  when  applying 
this  approach  to  tissues  with  a  markedly  high  macromolec¬ 
ular  content.  Additionally,  a  combination  of  either  mCRI  fit 
or  the  first-order  correction  with  the  recent  single-point 
CRI  technique  (12)  provides  a  promising  approach  for  the 
fast  unbiased  estimation  of  water  proton  Rt  values  within 
the  VFA  method  (35).  Finally,  it  should  be  pointed  out  that 
other  sources  of  bias  in  the  VFA  method  need  to  be  elimi¬ 
nated  in  conjunction  with  the  proposed  methodology,  simi¬ 
lar  to  the  procedures  described  in  this  study.  Practically, 
VFA  R1  mapping  should  be  performed  with  a  3D  sequence 
providing  a  uniform  slab  profile,  B1  correction,  and  appro¬ 
priate  spoiling  conditions  (30,31). 

The  described  analysis  and  correction  methods  have 
several  potential  limitations.  One  is  the  fundamental  li¬ 
mitation  of  the  two-pool  model  for  accurate  description 
of  multiple  exchanging  compartments  in  biological  tis¬ 
sues  (36,37).  Another  is  a  specific  assumption  about  the 
longitudinal  relaxation  rate  of  the  bound  pool  Rf .  Partic¬ 
ularly,  a  recent  study  (26)  demonstrated  that  widely 
adapted  approximations  =  1  s-1  (2,5)  or  R^  =  Rf 
(6,12)  may  considerably  underestimate  actual  values  of 
this  parameter.  Accordingly,  this  may  cause  an  addi¬ 
tional  bias  in  the  two-pool  model  parameters  obtained 
using  all  qMT  techniques  based  on  the  above  assump¬ 
tions  about  Rf.  However,  it  remains  unclear  at  this  point 
whether  a  faster  longitudinal  relaxation  rate  of  the  bound 
pool  protons  is  an  intrinsic  property  of  macromolecular 
protons  or  a  manifestation  of  more  complex  equilibria 
involving  water  proton  fractions  with  different  mobilities. 
The  later  hypothesis  is  logical  in  view  of  a  shorter  T1 
found  for  the  short-T2  water  fraction  commonly  termed 
“myelin  water”  (36,38-41).  However,  another  limitation  is 
the  assumption  about  a  negligible  effect  of  the  excitation 
pulse  on  the  bound  pool.  This  approximation  is  common 
for  off-resonance  pulsed  quantitative  MT  methods 
employing  low-angle  excitation  pulses  (5-7,10,12),  which 
causes  negligible  on-resonance  saturation  of  the  bound 
pool  magnetization.  However,  as  was  shown  earlier  (25), 
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the  effect  of  on-resonance  saturation  should  be  taken  into 
account  for  sequences  with  very  short  TR  and  high  flip 
angles,  the  conditions  typical  for  the  balanced  steady  state 
precession  imaging.  In  the  two-pool  formalism  with 
pulsed  excitation  (Eq.  [1]),  on-resonance  saturation  of  the 
bound  pool  can  be  taken  into  account  by  assigning  an 
actual  saturation  factor  Sb  <  1  in  the  diagonal  matrix  C. 
Corresponding  Sb  values  can  be  estimated  from  the  prop¬ 
erties  of  the  excitation  pulse  (duration  and  flip  angle)  and 
a  spectral  line  shape  of  the  bound  pool  (11,22).  Based  on 
this  approach,  we  found  for  a  particular  pulse  used  in  our 
sequence  that  on-resonance  saturation  of  the  bound  pool 
results  in  an  almost  negligible  effect  on  the  accuracy  of 
the  two-pool  model  parameters  determined  by  mCRI 
(<1%  relative  error  for  comparison  between  fits  with 
actual  Sb  vs.  Sb  =1).  However,  this  aspect  of  the  tech¬ 
nique  could  be  important  in  the  context  of  its  standardiza¬ 
tion  for  future  multicenter  clinical  studies  since  variations 
in  the  shape,  duration,  and  flip  angle  of  the  excitation 
pulse  may  introduce  a  small  but  non-negligible  bias  de¬ 
pendent  on  a  particular  sequence  implementation.  Simi¬ 
larly,  the  approximation  of  the  effect  of  the  excitation 
pulse  on  the  free  pool  by  the  term  cos  a  in  the  matrix  C 
(Eq.  [1])  may  not  be  absolutely  correct  due  to  relaxation 
during  RF  pulses  of  finite  duration.  Such  an  effect  has 
been  recently  described  (42),  and  its  role  for  both  VFA  R1 
and  CRI  parameter  measurements  remains  to  be  investi¬ 
gated.  In  summary,  while  this  study  relies  on  the  com¬ 
monly  accepted  in  the  quantitative  MT  methodology 
approximations,  future  refinements  of  CRI  accuracy  seem 
to  be  possible  based  on  accommodation  of  multicompart¬ 
ment  magnetization  exchange  schemes  and  a  more  rigor¬ 
ous  analysis  of  magnetization  dynamics  during  RF  pulses. 

CONCLUSIONS 

This  study  demonstrates  on  the  example  of  CRI  that  the 
accuracy  of  quantitative  MT  imaging  can  be  considerably 
improved  by  taking  into  account  the  contribution  of  MT- 
induced  biexponential  longitudinal  relaxation  into  VFA 
R1  measurements.  We  have  developed  the  two  technical 
approaches  allowing  correction  of  biases  in  CRI  parameter 
estimates  caused  by  the  MT-induced  biexponential  behav¬ 
ior  of  longitudinal  relaxation  of  water  spins  based  on  ei¬ 
ther  the  global  fit  of  both  VFA  and  MT  SPGR  data  to  the 
two-pool  signal  model  or  the  simple  analytical  recalcula¬ 
tion  formulas  for  key  cross-relaxation  parameters. 
Although  the  choice  between  these  techniques  depends 
on  a  tradeoff  between  the  correction  accuracy  and  image 
processing  speed,  both  approaches  do  not  require  extra 
data  compared  to  the  original  CRI  method. 
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Fig  1:  Diagram  of  BS-AFI  pulse  sequence,  a  is  a  slab-selective 
excitation  pulse,  while  ±co  is  a  nonselective  Fermi  pulse  applied 
at  ±co  Hz  with  respect  to  center  frequency.  Additional  Fermi 
pulses  continue  to  play  at  regular  intervals  during  the  second 
TR  (figure  compressed  to  save  space). 
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Fig  2:  Comparison  of  standard  AFI  (a)  &  BS  (d)  methods  to  the 
combined  BS-AFI  approach  (b,e),  normalized  by  nominal  flip 
angle.  Significant  differences  between  FA  (a,b)  and  B1  (d,e) 
are  observed,  primarily  due  to  slab  profile  (vertical  direction). 
Absolute  error  (c,f)  is  quite  uniform,  despite  large  off-resonance 
within  slice  (-200  Hz,  not  shown). 
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Introduction:  The  accuracy  of  many  quantitative  MRI  (qMRI)  methods 
depends  critically  on  accurate  knowledge  of  radiofrequency  transmit  field 
amplitude  (Bj)  or  flip  angle  (FA).  While  the  terms  B,  and  FA  are  sometimes  used 
interchangeably,  FA  depends  not  only  on  Bj  but  also  additional  factors  such  as 
RF  pulse  envelope,  slab  profile,  tissue  relaxation  parameters,  and  B0  effects  that 
can  lead  to  substantial  errors  when  correcting  qMRI  methods  [1].  In  general,  Bj 
correction  is  most  appropriate  for  nonselective  pulses  (e.g.  off-resonance  pulse  in 
magnetization  transfer  (MT)  or  a  hard  180°  pulse  in  inversion  recovery)  while  FA  correction  is 
appropriate  for  selective  excitation  pulses  (e.g.  variable  flip  angle  77  mapping).  Many  methods  such  as 
quantitative  MT  [2,3]  and  Look-Locker  77  mapping  utilize  both  types  of  pulses  in  the  same  experiment, 
however  to  date  no  consideration  has  been  made  of  separate  B}  and  FA  effects.  Therefore,  the 
measurement  and  correction  of  both  parameters  may  increase  the  overall  accuracy  of  such  methods. 

Bloch-Siegert  Bi  mapping  (BS)  is  a  recently  proposed  phase-based  method  to  rapidly  map  7?y  with  an 
off-resonant  nonselective  pulse  [4],  while  Actual  Flip-angle  Imaging  (AFI)  is  a  method  to  map  FA 
using  the  ratio  of  signals  from  two  interleaved  steady-state  TRs  [5].  We  propose  a  novel  combination  of 
these  two  sequences,  BS-AFI,  to  map  FA  using  image  magnitude  and  encode  Bi  into  the  image  phase 
with  zero  time  penalty  over  AFI.  Improved  MT  accuracy  is  demonstrated  in  phantom  validations. 

Methods:  The  BS-AFI  pulse  sequence  was  implemented  on  a  GE  MR750  (fig.  1).  A  standard  AFI 
pulse  sequence  was  implemented  by  the  addition  of  a  second  TR2  (=5*TRi)  to  an  SPGR  sequence;  BS 
mapping  was  then  added  through  the  introduction  of  an  8ms  Fermi-shaped  off-resonance  pulse  just 
after  the  excitation  pulse  and  just  before  phase  encoding.  As  BS  requires  the  phase  differences  of  two 
±co  off-resonance  pulses,  a  pulse  at  +4kHz  is  played  during  TRi  and  -4kHz  during  TR2.  Four  additional 
alternating  pulses  are  played  during  TR2  to  maintain  a  pulsed  steady-state  saturation.  A  potential  issue 
of  this  implementation  is  if  MT  saturation  effects  are  different  between  TRi  and  TR2  due  to  variation  is 
center  frequency  with  B0.  However,  since  the  TR  of  the  AFI  sequence  is  shorter  than  the  timescale  of 
MT  exchange,  we  hypothesize  that  these  alternating  pulses  will  create  an  averaging  effect  and  generate 
the  same  MT  contrast  in  both  images,  eliminating  any  confounds  to  AFI  quantification.  The  Br 
dependant  phase  precession  may  also  interfere  with  RF  spoiling,  however  it  was  shown  that  the  use  of 
large  spoiler  gradients  eliminates  sensitivity  to  RF  spoiling  effects  [6]. 

A  cylindrical  gelatin  phantom  was  scanned  at  64x48x52  matrix  and  1.53  mm3  isotropic  resolution  in  the 
coronal  plane,  then  re-formatted  to  axial  to  visualize  slab  profile  (vertical  direction).  BS-AFI  and 
standard  AFI  scans  were  acquired  with  TRi/TR2  =  36/1 80ms  and  a  =  55°.  Spoiled  gradient  echo 
(SPGR)  scans  were  acquired  with  TR  =  36ms  and  a  =  [5  10  20  30  40]°.  Ten  MT-weighted  SPGR  were 
acquired  with  TR  =  36ms,  a  =  10°,  and  offset  frequency  A=[l  2. 5  49  13]  kHz  for  MT  flip  angles  500° 
and  920°.  The  SPGR  and  MT  data  were  then  fitted  to  a  five-parameter  model  of  MT  (bound  pool 
fraction,  exchange  rate,  T2  of  bound  pools,  PD,  and  77  of  free  pool)  [7].  The  qMT  data  were  corrected 
in  three  ways:  using  only  the  FA  map  to  correct  both  excitation  and  MT  pulse  flip  angles  (FA  only), 
using  only  the  B2  map  to  correct  both  (Bj  only),  and  using  the  BS-AFI  FA  to  correct  the  excitation  flip 
angle  and  Bj  to  correct  MT  flip  angle.  B0  correction  was  done  using  IDEAL  [8]. 

Results:  Mean  error  for  the  BS-AFI  method,  compared  to  separate  AFI  and  Bloch-Siegert  mapping, 
was  3.0  ±  2.2%  in  B1  and  2.0  ±  1.1%  in  FA  (fig.  2).  Spatially  varying  differences  between  techniques 
were  not  observed,  even  in  the  presence  of  B0  inhomogeneity  (0-200  Hz  within  slice,  as  measured  with 
IDEAL).  SNR  of  the  B}  map  was  also  improved  (variance  oStandard  =  0.0048  G,  c>bs-afi  =  0.0037  G; 
small  ROI  in  center  of  image)  due  to  increased  signal  from  a  longer  AFI  TR2.  qMT  parameter  maps 
using  only  FA  or  Bj  show  distinct  errors  due  to  slab  profile  (fig.  3,  see  caption).  BS-AFI  maps  show 
uniform  measures  along  the  slab  direction  (fig.  4). 

Discussion  &  Conclusions:  Simultaneous  knowledge  of  Bj  and  FA  may  be  useful  for  a  wide  range 
of  qMRI  applications  that  employ  both  selective  and  non-selective  RF  pulses,  including  quantitative 
MTI,  Look-Locker  77  mapping,  and  the  correction  of  stimulated  echoes  in  a  multi-exponential  T2 
CPMG  experiment  [9].  Although  errors  in  FA  due  to  slab  profile  can  be  minimized  in  other  ways,  such 
as  the  use  of  a  longer  optimized  RF  pulse  or  nonselective  3D  imaging,  such  approaches  are  not  always 
feasible.  Nonselective  imaging  is  only  acceptable  in  applications  where  the  object  is  sufficiently  small 
to  be  covered  by  two  phase  encoding  directions  in  a  reasonable  scan  time,  and  may  lead  to  additional 
errors  due  to  increased  sensitivity  to  B0  variations.  Pulses  that  are  optimized  for  slab  selection  may  still 
have  residual  errors  in  a  significant  number  of  slices,  as  we  have  observed  errors  greater  than  5%  in 
quantitative  maps  occur  in  1/3  of  our  slices.  These  approaches  also  do  not  consider  variations  in  FA  due 
to  B0  or  the  relaxation  properties  of  tissue  (e.g.  77  relaxation  in  ultrashort  TE  sequences). 

In  contrast,  the  BS-AFI  method  enables  rapid  and  accurate  qMRI  with  sequences  employing  any 
combination  of  selective  or  non-selective  RF  pulse  types.  The  key  feature  of  the  pulse  sequence  is 
encoding  the  two  alternating  off-resonance  shifts  into  two  AFI  TRs.  This  may  potentially  influence 
contrast  of  AFI  images  differently  if  TR  is  long,  however  no  variations  were  observed  in  our  experiment  using  typical  AFI  parameters.  [3,6].  This  averaging  effect  may 
be  beneficial  not  only  in  the  context  of  BS-AFI,  but  also  ^^-insensitive  qMT  or  MT  ratio  measurements,  and  will  be  the  subject  of  future  work.  In  conclusion,  we 
present  a  method  to  simultaneously  map  B1  and  FA  that  is  robust  to  B0  and  MT  errors  and  incurs  no  scan  time  penalty  over  the  traditional  AFI  method.  In  context  of 
MT,  separate  measurement  of  FA  for  excitation  flip  angle  correction  and  BI  for  MT  power  correction  led  to  a  significant  increase  in  the  accuracy  of  bound  pool 
fraction.  In  addition  to  enlongated  scan  times,  separate  acquisition  of  such  maps  may  result  in  reduced  consistantcy  of  Bj  and  FA  mapping  due  to  patient  motion. 

References:  [1]  Wang  et  al.  MRM  2006;56:463  [2]  Sled  JG  et  al.  MRM  2001;46:923  [3]  Yamykh  VL  MRM  2002;47:929  [4]  Sacolick,  LI  et  al.  MRM  2010;63:1315 
[5]  Yamykh  VL  MRM  2007;57:192.  [6]  Yamykh  VL  MRM,  2010;63:1610  [7]  Mossahebi  P  et  al.  WMSG  2011  pl3.  [8]  Yu,  H.  JMRI  2007;26:1 153  [9]  Prasloski  T.  et 
al.  MRM  201 1  (in  press).  Acknowledgements:  We  acknowledge  the  financial  support  of  the  NIH  (NINDS  R01NS065034). 


Fig  3:  R1  (=1/77)  (a-c)  and  bound  pool  fraction  f  (d-f)  maps 
using  three  correction  approaches.  Note  that  R1  mapping  using 
B-i  correction  (b)  or  correction  of  MT  pulse  with  FA  (d)  both 
result  in  large  residual  errors  due  to  slab  profile.  Similar  errors 
were  observed  in  PD  and  exchange  rate  (not  shown). 


Fig  4:  Plot  of  R1  and  f  along  slab  profile  (dash  line  in  fig.  3a). 
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INTRODUCTION:  Over  past  years,  several  methods  have  been  proposed  for  quantitative  mapping  of  molecular  parameters  describing  the  magnetization  transfer 
(MT)  effect  within  the  two-pool  model  (bound  pool  fraction  f  cross-relaxation  rate  k,  and  transverse  relaxation  times  of  the  bound  pool  T2B  and  free  pool  T2F)  from  a 
series  of  images  with  variable  off-resonance  saturation  [1-4].  Resulting  parametric  maps  are  expected  to  provide  new  biomarkers  in  various  neurological  conditions. 
Particularly  encouraging  results  were  reported  for  the  bound  pool  fraction,  which  provides  strong  associations  with  myelin  content  in  neural  tissues  [5,6].  The  common 
feature  of  the  above  methods  is  the  need  for  complementary  T1  mapping,  which  allows  decoupling  of  the  two-pool  model  parameters  from  the  longitudinal  relaxation 
rate  Rl=l/Tl.  One  of  such  methods,  cross-relaxation  imaging  (CRI)  [4]  allows  3D  acquisition  with  excellent  speed  and  resolution  based  on  combination  of  a  limited 
number  of  MT-weighted  images  and  variable  flip  angle  (VFA)  T1  maps  obtained  using  a  fast  spoiled  gradient-echo  (SPGR)  sequence.  However,  recent  studies  showed 
that  the  on-resonance  MT  effect  may  cause  a  significant  bias  in  the  steady-state  MR  signal  [7,8],  if  the  model  does  not  account  for  bi-exponential  behavior  of 
longitudinal  relaxation  occurring  due  to  MT.  This  bias  unavoidably  introduces  errors  in  VFA  R1  measurements  [8],  which  may  propagate  into  the  two-pool  model 
parameters  measured  by  CRI.  In  this  study,  we  theoretically  and  experimentally  characterize  systematic  errors  in  CRI  caused  by  MT-related  bias  in  VFA  R1  mapping 
and  propose  a  modified  processing  algorithm,  which  corrects  for  such  errors  and  yields  accurate  quantitative  MT  parametric  maps.  Modified  Original 

THEORY:  In  the  standard  CRI  approach,  reconstruction  of  parametric  maps  is  performed  in  two  steps  [4].  During  the  first  step,  R1 
and  proton  density  {PD)  maps  are  calculated  from  VFA  data  by  fitting  the  Ernst  equation.  This  approach  ignores  coupling  between 
R1  and  the  rest  of  two-pool  model  parameters.  During  the  second  step,  the  matrix  model  of  pulsed  MT  [4]  is  fitted  to  MT-weighted 
data,  while  R1  is  supplied  as  an  external  parameter.  Prior  to  fit,  MT  data  are  normalized  to  the  signal  intensity  obtained  without 
saturation  to  exclude  PD.  Depending  on  the  sampling  scheme  and  constraints  imposed  on  certain  two-pool  model  parameters,  either 
all  four  parameters  (f,  k,  T2B,  T2F)  or  their  subsets  can  be  fitted  [4,5,9].  As  an  alternative  to  this  approach,  we  modified  the  CRI 
model  to  simultaneously  perform  global  fit  of  VFA  and  MT  data,  which  yields  all  parameters  (PD,  Rl,  f  k,  T? ,  T2F)  with  full 
consideration  of  coupling  between  them.  Theoretically,  MR-related  errors  in  Rl  measured  by  the  VFA  method  can  be  caused  by  the 
two  mechanisms  unaccounted  within  the  standard  single-pool  signal  model.  These  are  on-resonance  saturation  of  the  bound  pool  by 
the  excitation  pulse  and  bi-exponential  relaxation  during  TR,  where  a  contribution  of  the  fast-relaxing  component  becomes  more 
significant  at  short  TR.  To  estimate  relative  roles  of  these  factors,  the  model  was  further  modified  to  include  an  additional  term  to 
account  for  the  effect  of  excitation  pulse  on  semisolid  pool  [7]  calculated  by  Bloch  equation  simulations. 

METHODS:  Volunteer  data  were  obtained  on  a  3.0T  GE  MR750  (Waukesha,  WI)  scanner.  To  correct  for  B1  and  BO 
inhomogeneity,  flip  angle  (Bl)  and  field  (BO)  maps  were  measured  by  optimized  AFI  method  [10]  and  IDEAL  [11],  respectively. 

3D  VFA  data  were  acquired  with  FAs  a=[5, 10,20,30]°.  Ten  MT  SPGR  datasets  were  acquired  (a=10°,  A  =0.8,  2.5,5,9,13  kHz, 
oimt=[500,1  100]°,  18ms  Fermi  pulse).  Additional  dataset  without  saturation  was  acquired  at  offset  frequency  200kHz  to  normalize 
MT  data  for  the  original  CRI  approach.  All  data  were  acquired  with  TR/TE=40/2.0ms,  240x  180><80mm  FOV,  128x96x40  matrix. 

The  standard  CRI  processing  workflow  was  implemented  according  to  [4,9].  In  modified  CRI,  all  parametric  maps  were  generated 
by  fitting  the  SPGR  and  MT  data  simultaneously  to  the  general  signal  equation  based  on  the  two-pool  matrix  steady-state  MT  model 
[2,4]  using  in-house-written  C  and  MATLAB  (MathWorks,  Natick,  MA)  software  utilizing  a  trust-region-reflective  algorithm  for 
nonlinear  least  squares  fitting  of  each  voxel. 

RESULTS:  Synthetic  MT  and  SPGR  data  for  white  matter  (WM)  and  gray  matter  (GM)  ROIs  were  generated  based  on  the 
parameters  given  in  [9].  Table  1  shows  that  on-resonance  MT  effect  introduces  significant  bias  to  all  quantitative  MT  parameters 
except  T2b  when  estimated  with  original  CRI  method.  The  dominant  cause  of  errors  was  found  to  originate  from  effective  shortening 
of  T1  in  the  two-pool  system  due  to  bi-exponential  relaxation,  whereas  saturation  of  the  bound  pool  by  the  excitation  pulse  produced 
an  almost  negligible  effect  (~1%).  Qualitative  (Table  2)  and  quantitative  in  vivo  comparisons  (Fig.  1)  of  original  and  modified  CRI 
methods  are  consistent  with  simulations  and  demonstrate  compatible  effect  sizes  in  all  parameters. 

DISCUSSION:  We  showed  that  errors  caused  by  unaccounted  effect  of  MT  on  VFA  Rl  measurements  result  in  significant  bias  in 
estimates  of  CRI  parameters.  Generally,  this  bias  depends  on  the  values  of  the  two-pool  model  parameters,  which  makes  it 
inconsistent  for  different  tissues  types  and  in  different  pathological  conditions.  Approximately,  this  type  of  errors  can  be 
characterized  as  overestimation  of  Rl,  f  and  k,  and  underestimation  of  T2F  with  a  relative  effect  size  on  the  order  of  /%.  Our  results 
show  that  accuracy  of  the  original  CRI  approach  is  noticeably  improved  with  the  proposed  combined  data  fit  to  the  modified  CRI  Figure  1.  Parametric  MT  maps 

model,  which  accounts  for  MT  effects  on  the  apparent  Rl.  In  brain  imaging,  the  method  also  corrected  VFA  Rl  values  by  estimated  using  modified  (left) 

approximately  15%  in  WM  and  6%  in  GM,  which  agrees  well  with  the  bias  predicted  by  [8].  Finally,  the  modified  CRI  fit  does  not  and  original  (right)  CRI. 
require  acquisition  of  normalizing 

CONCLUSIONS:  This  study  demonstrates  that  separate  treatment  of  VFA  and  MT  data 
in  the  CRI  method  causes  non-negligible  systematic  errors  in  both  Rl  and  cross¬ 
relaxation  parameters.  Our  modified  CRI  data  processing  approach  effectively  corrects 
these  errors  and  does  not  require  any  additional  measurements,  thus  maintaining  time- 
efficiency  of  the  original  CRI  technique. 


Modified  CRI 

Original  CRI 

Rl  (s1) 

f(%)  k(s') 

T2b  (jus) 

T2f  (ms) 

Rl  (s')  f(%) 

k  (s'1) 

T2b  (jus) 

T2f  (ms) 

WM  (Corpus  Callosum) 

0.977±0.032 

14.00±0.46  2.635±0.218 

10.31+0.32 

25.06±2.23 

1.104±0.038  16.62+0.78 

2.815±0.242 

10. 15+0.32 

21.77+1.95 

GM  (Putamen) 

0.726±0.039 

7.21±0.39  1.359±0.094 

10.40±0.31 

33.06±2.89 

0.774±0.044  8.12±0.49 

1.363+0.098 

10.24±0.30 

30.95±2.70 

Table  2.  In  vivo  measurements  of  CRI  parameters  using  modified  and  original  CRI. 
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Rl 

/ 

k 

T2b 

T2f 

WM  (Corpus  Callosum) 

15.27% 

11.24% 

13.97% 

2.76e-4% 

-11.49% 

GM  (Putamen) 

8.16% 

5.65% 

6.49% 

2.63e-4% 

-5.64% 

Table  1.  Errors  in  two-pool  MT  model  parameters  in  original  CRI 
fit  due  to  unaccounted  on-resonance  MT  effects  (simulations). 
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INTRODUCTION:  Quantification  of  two-pool  MT  (qMT)  model  parameters  (longitudinal  relaxation  rate  Rl,  bound  pool  fraction  f  cross  relaxation  rate  k,  transverse 
relaxation  times  of  the  bound  pool  T2B,  and  free  pool  Tf)  allows  inferring  useful  information  about  MRI-invisible  macromolecular  tissue  content  such  as  myelin  [1,2]. 
Many  qMT  methods  such  as  cross-relaxation  imaging  (CRI)  [3]  use  efficient  pulsed  MT  in  steady  state,  however  data  acquisition  time  may  be  too  long  for  clinical 
applications.  One  CRI  acceleration  approach  is  to  exclude  Tf  from  estimation  by  sampling  high  offset  frequencies  (A>2.5  kHz)  where  direct  saturation  of  free  water  is 
minimized.  Going  from  1.5T,  the  system  used  in  early  qMT  applications  [1-3],  to  3T  may  improve  CRI  by  utilizing  the  field-dependent  SNR  increase  [4].  However, 
translation  to  3T  faces  an  increased  specific  absorption  rate  (SAR);  the  related  decrease  in  highest  attainable  MT  flip  angle  may  lead  to  SNR  loss  in  parametric  maps. 
This  work  studies  challenges  of  CRI  at  3T  and  proposes  several  optimizations  at  both  1.5  and  3T  field  strengths. 

METHODS  AND  RESULTS:  Effect  of  Experiment  Desisn:  Two  families  of  experiment 
designs  were  studied:  1)  Standard  CRI  with  higher  A=2. 5,  5,  9,  13  kHz  (only  Rl,  f  k,  T2B  are 
estimated).  2)  Extended  design  with  additional  optimized  low  A  =0.8  kHz  {Tf  is  also 
estimated).  Data  were  synthesized  for  white  matter  for  18ms  Fermi  MT  pulse  with  flip  angles 
500°  (1.5T  and  3T)  and  1100°(3T)  /1600°(1-5T)  (maximal  allowed  values  on  our  GE  Signa 
1.5T  and  MR750  3T  for  the  given  patient  weight  and  coil  type).  Monte-Carlo  simulations  were 
performed  to  find  noise  level  o  in  qMT  estimates  (Fig.  1).  Efficiency  was  defined  as  1/a  and 
did  not  include  field  dependent  SNR  gain.  Using  low  A  significantly  improved  estimation  of 
CRI  measures  (especially  k)  despite  dimensionality  increase  {Tf  estimation).  Effect  of  MT 
Pulse  Shape:  The  effect  of  using  MT  pulses  with  low  (Gaussian)  and  high  (Fermi)  SAR- 
efficiency  [6]  on  the  qMT  mapping  performance  were  examined  for  Design  #2.  The  pulses 
had  same  duration  and  peak  SAR.  Fermi  pulse  was  more  SNR-efficient  for  parameter 
estimation  than  Gaussian  pulse  with  most  gain  observed  for  k  (Table  1).  Effect  of  MT  Pulse 
Bandwidth:  Higher  BW  (shorter)  Fermi  pulse  may  cause  large  errors  in  qMT  maps  using  low 
offset  frequencies  (Table  2).  The  standard  approach  is  to  use  Lorentzian  lineshape  to  model  the 
effect  of  MT  pulse  on  free  pool.  The  approach  fails  to  adequately  describe  the  behavior  of  free 
pool  in  the  range<lkHz  (Fig.  2c).  The  deviation  is  due  to  Rabi  oscillations  arising  in  near¬ 
resonance  excitation  regime,  the  effect  highly  dependent  on  local  B1  distribution  (Fig.  2a)  not 
accounted  for  in  realistic  qMT  models  [5].  The  MT  signal  is  much  more  sensitive  to  changes 
in  BO  and  B1  for  shorter  pulse  (Figs.  2a, b).  The  effect  is  hard  to  account  for  in  realistic  qMT 
models  [5]  because  it  is  highly  dependent  on  local  distribution  of  transmit  field  B1  (Fig.  2). 

The  accuracy  loss  for  20  ms  pulse  (narrower  BW)  was  much  less  for  most  measures  (Rl,  f  k, 

T2b )  while  still  significant  for  Tf  (Table  2).  Phantom  Results:  MT  weighted  data  of  FBRIN 
phantom  [7]  were  acquired  at  offset  frequencies  0.6  and  20kHz  with  8ms  Fermi  MT  pulse  with 
flip  angle  750°.  Range  of  offset  frequency  in  the  phantom  was  provided  using  shim  gradients 
and  measured  using  IDEAL  [8].  Flip  angle  map  was  measured  by  optimized  AFI  method  [9]. 

The  3D  images  of  MT  ratio  (SMT(0.6kHz)/SMT(20kHz))  were  plotted  against  local  flip  angle 
(~B1)  and  B0  values  (Fig.  2d).  The  observed  oscillations  were  consistent  with  simulations  for 
free  water  magnetization  (Fig.  2b)  (full  quantitative  comparison  between  theoretical 
predication  phantom  data  was  not  made  because  of  additional  MT  effect  in  the  agar  phantom 
and  uncertainties  in  B1  field  determination  from  actual  flip  angle  values). 

DISCUSSION:  We  studied  several  ways  to  optimize  efficiency  of  CRI  at  both  1.5  and  3T.  We 
showed  (for  3T  case)  that  the  CRI  efficiency  may  be  improved  by  including  extra  sample(s) 
from  low  A  range  (<1  kHz).  Further,  we  demonstrated  that  SAR-efficient  pulses  may  provide 
additional  large  efficiency  gains  in  CRI.  At  the  same  time,  high  BW  SAR-efficient  pulses  may 
lead  to  a  significant  error  in  qMT  measures  for  the  low  offset  frequency  designs.  One  simple 
solution  is  to  resort  to  longer  SAR-efficient  pulses.  We  found  the  reasonable  tradeoff  to  be  lb- 
20  ms  Fermi  pulses.  We  are  currently  working  on  MT  pulse  BW  correction  for  more  efficient 
imaging  with  shorter  pulses.  We  conclude  that  CRI  and  other  related  qMTI  methods  may  be 
optimized  as  proposed  here  to  significantly  improve  their  accuracy  and  efficiency. 
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Fig.  1:  Relative  efficiency  of  CRI  parameter  mapping  at  3T  using 
regular  (#1)  and  extended  (#2)  designs  (normalized  to  1.5T 
values). 
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Fig.  2:  Normalized  MT  signal  vs.  offset  frequency  and  MT  flip 
angles  for  8ms  (a)  and  20ms  (b)  MT  pulse  widths  shows 
significant  Rabi  oscillations  for  higher  BW  pulse.  The  Rabi 
oscillations  make  the  standard  Lorentian  lineshape  approximation 
inadequate  to  model  effect  of  higher  BW  pulse  on  free  pool  for 
low  offset  frequencies  (c).  MT  signal  from  FBRIN  phantom  for 
the  range  of  A  and  ocMt. 


Rl 

f 

k 

T  B 

1  2 

t2f 

1.0413 

1.5820 

2.9149 

1.7558 

1.9361 

Table  1 :  Efficiency  of  Fermi  relative  to  that  of 
GaUSSian  MT  pulses  (Ggauss/tffermi)- 


Rl 

f 

k 

/r1 

CD 

8ms  Fermi 

-11.33 

40.7 

-48.7 

18.23 

20ms  Fermi 

-0.01 

-0.12 

-3.24 

-0.32 

Table  2:  Percent  of  error  in  CRI  parameters  for 
different  MT  pulse  widths  (Design  #2). 
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Figure  1.  a:  Cortical  thickness  estimated  from  IR-SPGR  and  MPF  (both  1  mmA3).  IR- 
SPGR  underestimates  the  thickness  in  areas  of  reduced  WM/GM  contrast  arising  due  to 
Bi  field  variations  typical  for  high  field  MRI.  Note  increased  spatial  uniformity  of 
cortical  thickness  measured  with  the  quantitative  myelin-sensitive  MPF.  The  Bl  effect 
is  illustrated  in  (b,  c).  Reduced  IR-SPGR  WM/GM  contrast  in  the  right  hemisphere  (b) 
corresponds  to  the  areas  of  elevated  Bi  (c). 
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Figure  2.  Improved  visualization  of  deep  GM  structures  on  high  resolution,  myelin- 
sensitive  MPF  compared  to  the  resolution-matched  quantitative  Ri  and  inhomogeneity 
corrected  anatomical  IR-SPGR  maps  (arrows). 


Introduction:  Two-pool  modeling  of  magnetization  transfer  (MT)  effects  yields  a  unique  set  of  measures  sensitive  to  different  tissue  composition  properties  [1-4].  The  key 
parameter  of  interest  in  this  model  is  a  macromolecular  pool  fraction  (MPF),  a  potential  biomarker  of  myelin  in  neural  tissues  [5-7].  Particularly  promising  findings  were  recently 
reported  on  the  utility  of  MT-based  measures  including  MPF  for  assessment  of  white  and  gray  matter  (WM/GM)  myelin  status  in  multiple  sclerosis  (MS)  [7-8].  In  particular,  MPF 
in  GM  was  found  to  be  the  strongest  independent  predictor  of  disability  in  MS  patients  [7,9].  Given  that  cortical  demyelination  has  been  recognized  as  an  important  pathological 
substrate  in  MS  [10,1 1],  the  availability  of  such  a  myelin  biomarker  may  greatly  facilitate  our  understanding  of  GM  involvement  in  MS  disease  course.  However,  traditional  MPF 
mapping  within  clinical  scan  time  is  a  relatively  low-resolution  methodology  not  fully  compatible  with  current  standards  for  imaging  of  cortical  stmctures  (~1  mm  isotropic)  or 
achieving  higher  fidelity  in  detection  of  demyelinating  lesions  in  WM.  Longer  scan  times  needed  to  attain  the  resolution  makes  MPF  mapping  susceptible  to  small  patient  motion. 
The  purpose  of  this  work  was  to  develop  a  clinically  feasible,  high  resolution,  motion  corrected  MPF  mapping  protocol  compatible  with  clinical  standards  for  GM  assessment. 
Methods:  Three  Point  MPF  Manning:  Fast  cross-relaxation  imaging  (CRI)  [13]  was 
proposed  as  an  efficient  strategy  for  MPF  mapping  with  reduced  number  of  acquisition  points 
including  spoiled  gradient  echo  (SPGR)  with  (MTON)  and  without  (MT0ff)  MT  weightings, 
and  two  variable  flip  angle  (VFA)  SPGR  datasets  for  pre-requisite  Ri  mapping.  This  method 
imposes  realistic  constraints  onto  the  subset  of  two-pool  MT  model  parameters  (exchange  rate 
R  and  T2  of  bound  protons)  and  prescribes  an  optimized  experimental  design  to  provide  least 
sensitivity  to  these  constrained  quantities.  We  combined  the  constrained  approach  with  a 
modified  CRI  (mCRI)  fit  [13]  to  remove  the  bias  in  estimation  due  to  on-resonance  MT 
effects  and  also  further  reduce  the  number  of  measurements  by  25%  (from  four  to  three).  Our 
approach  (3p-mCRI)  performs  simultaneous  fitting  of  two  VFA  and  one  MT  SPGR  datasets 
to  the  unified  signal  model,  which  allowed  estimation  of  MPF,  Rj,  and  proton  density  terms 
without  the  MT0ff  dataset  [14].  3D  Undersamvled  Motion-Corrected  Radial  Acquisition:  To 
study  feasibility  of  high  resolution,  motion-corrected  MPF  mapping  in  clinically  feasible 
time,  the  method  was  implemented  with  3D  radial  undersampled  acquisition  trajectory.  3D 
radial  sampling  trajectory  provides  mild  undersampling  artifacts  allowing  accelerated  imaging 
without  noise-amplifying  matrix-inverse-based  parallel  MRI  techniques  (pMRI). 

Additionally,  a  notable  advantage  of  radials  is  the  ability  to  implement  interleaved  projection 
ordering  to  acquire  low-resolution  motion-tracking  navigators  from  the  oversampled  k-space 
center.  We  implemented  2D  bit  reversed  scheduling  of  projections  to  acquire  such  navigators 
to  enable  retrospective  motion  correction  as  described  in  [15].  As  non-uniform  &-space 
sampling  density  of  3D  radials  incurs  additional  scan  time  penalty  (a  factor  of  -1.57),  we 
designed  variable  density  readout  (VDS)  to  partially  compensate  for  this  efficiency  loss, 
which  theoretically  allows  improving  SNR  efficiency  by  -1.37  compared  to  uniform  readout 
trajectory  of  equivalent  length  [16].  Data  Acquisition:  Cartesian  Protocols:  SPGR: 

TR=23ms,  a=[4,24]°;  MT  SPGR:  TR=37ms,  a=12°,  A=4  kHz,  aMT=[850]°,  18ms  Fermi  pulse;  Common  parameters:  TE=3.7  ms,  2x  GRAPPA, 

8  channel  coil.  3D  Radials:  SPGR:  TR=18ms,  a=[4,22]°;  MT  SPGR:  TR=33ms,  a=13°,  A=4  kHz,  aMT=[750]°,  12ms  Fermi  pulse;  Common 
parameters:  22,000  projections  (undersampling  factor  R=3.45),  FOV=220,  TE=6.2  ms,  32  channel  coil.  All  protocols  included  AFI  flip  angle 
mapping  [17].  Total  scan  time  was:  15.5/29  min  (Cartesian  scans  with  1.25/1  mm  isotropic  resolution,  respectively),  and  28  min  for  radial 
scans.  Data  were  acquired  using  on  a  3.0T  GE  MR750  (Waukesha,  WI).  Processing  was  done  by  in-house-written  software  [18].  Cortical 
thickness  mapping  was  performed  with  FreeSurfer  software  [19]. 

Results:  One  healthy  volunteer  and  one  relapsing-remitting  (RR)  MS  patient  were  scanned  with  the  developed  protocols.  3p-mCRI  MPF 
mapping  was  first  validated  in  Cartesian  scans.  Figure  1  compares  cortical  thickness  measured  with  different  imaging  techniques.  3p-mCRI 
MPF  outperforms  the  standard  IR-SPGR  sequence  used  most  often  for  measurements  of  cortical  atrophy.  Figure  2  shows  improved 
visualization  of  deep  GM  stmctures  with  MPF  compared  to  IR-SPGR  and  Rj  maps  (all  1  mm  isotropic).  The  longitudinal  stability  of  3p-mCRI 
was  tested  in  a  RR  MS  patient  during  multiple  visits.  Figure  3  demonstrates  excellent  longitudinal  consistency  of  3p-mCRI  MPF  mapping. 

Figure  4a  shows  performance  of  3D  radial  implementation  of  3p-MRI  MPF  mapping.  The  protocol  provides  excellent  depiction  of  cortical  and 
deep  GM  structures.  Figure  4b  illustrates  its  motion  correction  capabilities  of  the  radial  protocol. 

Discussion:  We  presented  a  series  of  methodological  and 
technical  developments  aimed  towards  robust  methodology 
for  high  resolution  MPF  mapping.  Compared  to  early 
published  protocols  for  fast  MPF  mapping  [13,14],  our  3D 
radial  approach  decreases  voxel  size  by  a  factor  of  -7-11  at 
the  expense  of  -2.2-2. 5  times  scan  time  increase,  while 
maintaining  adequate  SNR  (Fig.  4)  and  clinically  acceptable 
scan  time  (<30  min).  The  imaging  protocol  derives  from 
several  sources,  including  minimalistic  3-point  constrained 
mCRI  approach,  SNR-optimized  3D  radials,  and  acquisition 
scheduling  for  motion  correction.  Acquisition  of 
undersampled  radials  with  32-channel  array  did  not  require 
specialized  pMRI  reconstruction  (well  known  property  of 
radial  sampling  at  low  accelerations)  and  hence  avoided  noise 
amplifications  common  in  matrix-in  verson  based  pMRI  techniques.  However,  the  dedicated  pMRI  reconstmction  may  be  required  when  imaging  is  performed  with  fewer  coil 
elements,  or  when  higher  acceleration  factors  are  targeted.  The  protocol  provides  motion-corrected  and  longitudinally  stable  MPF  (and  Rf)  maps,  which  may  establish  it  as  a  tool 
for  high  precision  morphometric  applications  and  GM/WM  segmentation  in  longitudinal  and  cross-sectional  studies  (Figs.  1,3).  Myelin-specific  contrast  obtained  with  MPF  may 
be  especially  important  to  study  diffuse  changes  in  both  GM  and  WM.  Reduced  partial  voluming  with  CSF  due  to  higher  resolution  may  be  important  to  provide  imaging 
assessment  of  subpial  GM  involvement  in  MS,  while  stronger  WM/GM  contrast  in  MPF  may  help  with  detection/classification  of  demyelinating  lesions  at  GM/WM  junction. 
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Figure  3.  Whole-brain  histograms 
of  MPF  obtained  by  3p-mCRI 
protocol  in  an  MS  subject  during 
5  monthly  visits. 


Figure  4.  Results  of  MPF  mapping  with  3D  undersampled  radials  at  lxlxl  mm  resolution,  a:  MPF  from  non-corrupted 
dataset  (two  axial  and  two  sagittal  slices  are  shown),  b:  MPF  from  motion-corrupted  dataset  before  and  after  motion 
correction  (two  coronal  slices  are  shown;  motion  caused  by  intentional  head  movement  of  the  healthy  volunteer).  Note 
significant  improvement  in  tissue  conspicuity  in  motion-corrected  MPF. 
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INTRODUCTION:  Variable  flip  angle  (VFA)  acquisition  using  spoiled  gradient  echo  (SPGR)  [1]  has  become  a  widely  accepted  approach  for  fast  high  resolution 
mapping  of  longitudinal  relaxation  rate  Rl=l/Tl.  Recent  studies  have  shown  that  magnetization  transfer  (MT)  between  free  and  bound  protons  due  to  the  application  of 
an  on-resonance  excitation  pulses  causes  a  significant  bias  in  steady  state  MR  signal  in  tissues  with  natural  abundance  of  macromolecular  content  [2].  Ou  et  al  predict 
up  to  14%  error  in  white  matter  (WM)  VFA  T1  values  [3],  an  effect  comparable  to  white  matter  T1  changes  observed  in  multiple  sclerosis  [4].  Explicit  consideration  of 
a  two-pool  MT  model  in  VFA  experiments  may  improve  the  accuracy  of  T1  estimation  [3].  However,  due  to  very  low  sensitivity  of  the  modified  VFA  signal  model  to 
parameters  of  the  two  pool  MT  model  such  as  bound  pool  fraction  (f),  cross-relaxation  rate  (k),  transverse  relaxation  times  of  the  bound  pool  T2B,  and  free  pool  T2F,  the 
feasibility  of  MT-corrected  in  vivo  T1  mapping  using  VFA  measurements  only  is  limited.  Cross-relaxation  imaging  (CRI)  provides  a  convenient  framework  to  model 
the  effect  of  both  off-resonance  MT  and  on-resonance  excitation  pulses  on  bound  and  free  proton  pools  [5].  It  was  recently  demonstrated  that  combination  of  VFA  and 
pulsed  MT  SPGR  measurements  within  modified  CRI  framework  may  be  used  to  remove  bias  due  to  on-resonance  MT  in  both  VFA  R1  values  and  two-pool  MT 
parameters  [6].  However,  this  approach  requires  augmenting  VFA  measurements  with  multiple  MT  acquisitions  to  allow  estimation  of  the  full  set  of  parameters.  Here, 
we  propose  a  method  to  obtain  accurate  T1  maps  of  clinical  quality  using  a  standard  VFA  datasets  augmented  with  a  single  off-resonance  MT  SPGR  scan  (VFA-MT). 

METHODS:  The  simultaneous  fit  of  SPGR  data  (VFA  part)  and  off-resonance  MT  SPGR  (MT  part)  data  to  the  modified  CRI  equation  accurately  yields  all  parameters 
of  the  model  (PD,Rl,fk,T2B)[ 6].  To  achieve  fast  77  mapping,  we  simplify  the  full  model  by  setting  T2B ,  T2fR f,  and  R=k(l-f)/f  to  constrained  values  as  proposed  in  [7] 
(11/rs,  0.022,  and  19s'1,  respectively),  which  leaves  only  three  free  parameters  {PD, Rif)  for  estimation.  These  parameters  may  be  in  principle  determined  using  only 
three  measurements  (two  VFA  and  one  MT).  The  accuracy  of  this  simplified  3 -parameters/3 -measurements  approach  (VFA-MT)  was  studied  in  simulations,  phantom 
imaging,  and  in  vivo.  Simulations:  We  run  simulations  to  identify  the  ranges  of  an  offset  frequency  and  an  MT  flip  angle  of  a  VFA-augmenting  MT  scan  for  which  the 
effect  of  such  constraints  on  accuracy  of  R1  and  /is  minimal.  Phantom  Studies:  Phantom  (25%  gelatin)  and  in  vivo  experiments  were  run  on  a  3.0T  GE  MR750 
(Waukesha,  WI),  and  flip  angle  and  field  maps  for  all  data  were  measured  by  an  optimized  AFI  method  [7]  and  IDEAL  [8],  respectively.  Phantom  protocol  included 
acquisition  of  3D  VFA  data  FAs=[7, 15,35,50]°  and  MT  SPGR  data  (FA=15°,  A  =2.5, 5, 9, 13kHz,  a  Mt=[850,1300]°,  18ms  Fermi  pulse,  120x120x88mm  FOV, 
128x128x44  matrix,  TR/TE=37/2.3ms).  Single-slice,  2D  inversion-prepared  spin-echo  (IR)  data  was  collected  to  determine  a  reference  T1  value  in  the  phantom 
(TR/TE=5 000/8 . 2ms ,  TI=0.05, 0.1, 0.2, 0.3, 0.5, 0.7, 0.9, 1.2, 1.6, 2s).  In  Vivo  Studies:  Volunteer  brain  data  included  3D  VFA  (FAs=[5, 10,20,30]°)  and  MT 
SPGR(FA=10°,  A  =2. 5, 5, 9, 13kHz,  aMT=[500,1100]°,  18ms  Fermi  pulse,  240x180x80mm  FOV,  128x96x40  matrix,  TR/TE=40/2.0ms).  The  reference  MT-corrected  PI 
maps  were  obtained  by  simultaneous  fit  of  all  VFA  and  MT  data  to  the  modified  CRI  model  [6].  All  processing  was  done  by  in-house  nonlinear  optimization  routines. 


RESULTS  AND  DISCUSSION:  Simulations:  Figure  1  illustrates  the  effect  of  two-pool  MT  model  parameters  on  observable  T1  values  when  calculated  using 
standard  VFA  method.  The  error  in  T1  increases  almost  linearly  with / for  typical  values  of  k  observed  in  neural  tissues  (>1)  and  hence  is  expected  to  be  most  biased  in 
macromolecular-rich  tissues  such  as  white  matter.  The  sensitivity  of  77  to  TB  was  negligible.  Figure  2  shows  the  effect  of  off-resonance  frequency  and  MT  flip  angle 
of  a  VFA-augmenting  MT  scan  on  the  average  loss  in  R1  and/for  the  proposed  simultaneous  constrained  estimation  of  these  parameters.  Similar  to  observations  in  [7], 
there  is  a  range  of  A  and  MT  flip  angle  with  low  sensitivity  to  the  fixed  parameters  (A=4-6kHz,  a  Mt=400-600°).  Therefore,  the  optimized  set  of  measurements  (SPGR 
a  =[5,30]°;  A=5kHz,  a  Mt=500°)  were  used  for  in  vivo  validation  of  the  constrained  model.  Phantom  Studies:  The  proposed  method  was  applied  to  correct  T1  mapping 
in  the  phantom  (f~ 4%).  Examination  of  mean  phantom  T1  values  obtained  with  full  CRI  (12  measurements),  proposed  constrained  VFA-MT  (3  measurements)  and 
VFA  (2  measurements)  confirmed  that  both  full  CRI  and  VFA-MT  provide  excellent  agreement  with  the  reference  IR  T1  measurement  (Fig.  3).  At  the  same  time,  VFA 
T1  was  biased  by  approximately /%,  which  is  consistent  with  simulations  in  Fig.  1.  In  Vivo  Studies:  Images  in  Fig.  4  reveal  that  there  is  a  visible  bias  in  T1  map 
obtained  with  regular  VFA  compared  to  T1  maps  obtained  with  on-resonance  MT  effect  correction.  The  MT-corrected  values  (77WM=1050+25,  77GM=1537+25)  agree 
well  with  literature  [10]  (77Wm~1060,  77gm~1630),  while  VFA  underestimated  them  (77WM=931+34,  T1GM= 1462+2 8).  The  whole  brain  T1  histograms  (Fig.  5)  show 
excellent  correlation  between  reference  full  CRI  values  and  proposed  VFA-MT  approach.  The  method  corrected  VFA  T1  values  by  11.5%  in  WM  and  6.2%  in  GM, 
which  agrees  well  with  the  bias  predicted  by  [3].  Analysis  of  ROI  77  measurements  (Table  2)  demonstrated  excellent  preservation  of  77  mapping  accuracy  in  both  WM 
and  GM  (error<l%)  when  moving  from  full  CRI  (12  measurements)  to  constrained  VFA-MT  (three  measurements)  estimation.  The  errors  in /maps  (not  shown  here) 
were  more  significant  but  always  <10%,  which  is  consistent  with  results  in  [7]. 

CONCLUSIONS:  Accurate  estimation  of  T1  values  from  VFA  measurements  requires  consideration  of  MT  effects,  which  may  introduce  T1  bias  roughly  proportional 
to  the  macromolecular  content.  The  fast  T1  correction  protocol  developed  in  this  work  (VFA-MT)  requires  a  single  MT  scan  in  addition  to  regular  VFA  measurements 
for  accurate  77  mapping.  The  presented  approach  is  promising  for  fast  high-resolution  whole  brain  T1  mapping  within  clinically  acceptable  scan  times. 
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Fig.  1  :  Error  of  VFA  T1  values  vs.  Fig.  2:  Dependence  of  error  in  estimation  of  R1=1/T1  (left)  and/(right)  due  to  constraining  vs. 

quantitative  MT  parameters /  and  k.  offset  frequency  and  MT  flip  angle  of  a  MT  SPGR  acquisition  of  proposed  3-point  protocol. 
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Fig.  3:  Phantom  T1  values 
calculated  by  different  methods. 


|3.5 
1 3.0 
5 

2.0 
1.5 

r 

1 0.0 

Fig.  4:  In  vivo  T1  maps  estimated  by  modified  full  (a) 
and  constrained  CRI  (b),  and  VFA  (c). 
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0.77% 

0.73% 

0.45% 

0.43% 

f 

3.05% 

5.89% 

4.22% 

9.73% 

Table  1:  Errors  of  constrained  method  in 
brain  structures  (CC:  Corpus  Callosum, 
IC:  Internal  Capsule,  PT:  Putamen,  CN: 
Caudate  Nucleus) 
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